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Abst ract 

A  development  and  analysis  of  a  single  controller, 
before  and  after  the  elimination  of  ''spillover”  terms,  is 
implemented  to  attempt  to  achieve  desired  response  charac¬ 
teristics  of  the  structure  under  evaluation.  Using  this 
derived  data  as  a  basis  for  comparison,  a  pair  of  decentral¬ 
ized  controllers  are  implemented  on  the  structure.  The 
characteristics  of  the  structural  response  is  dramatically 
improved  through  the  implementation  of  these  decentralized 
controllers.  Problems  encountered  with  the  implementation 
of  more  than  two  decentralized  controllers  are  investigated. 

The  structure  used  in  the  controller  evaluation  is  a 
lumped  mass  tetrahedron.  The  four  masses  of  this  model  are 
connected  through  isotropic  massless  rods  capable  of  support¬ 
ing  axial  loading  only  (no  bending) .  NASTRAN  is  used  to 
develop  a  normal  mode  approximation  of  the  structure,  while 
providing  mode  shape  and  frequencies  for  the  resultant 
twelve  mode  model.  Pointing  accuracy  of  the  apex  is  used  in 
determining  figures  of  merit  to  evaluate  the  effectiveness 
o f  the  control  applied.  Control  is  applied  through  each  of 
the  6  sensor/actuator  pairs  located  on  the  model . 

Controllers  are  developed  using  linear  optimal  tech¬ 
niques  which  produce  feedback  gains  proportional  to  the  state. 
The  state  is  represented  as  modal  amplitudes  and  velocities. 
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The  feedback  gains  are  established  via  steady  state  opt  imal 
rt  galat  or  theory.  The  system  response  is  evaluated  initially 
using  only  a  single  controller  on  an  eight  mode  truncated 
model  then  on  the  entire  twelve  mode  model.  A  comparison  is 
made  with  the  system  prior  to  the  eliminaf  ion  of  the  observa¬ 
tion  spillover  and  after  the  transformation  technique  is 
applied.  For  the  study,  four  modes  are  designated  as  control¬ 
led  and  four  as  suppressed.  The  remaining  modeled  modes  are 
designated  residual.  An  additional  controller  is  added  with  no 
addition  of  sensors  or  actuators.  While  the  response  of  the 
single  controller  system  is  unable  to  meet  the  design  criteria, 
the  addition  of  a  decentralized  controller  more  than  adequately 
achieves  the  desired  response. 


The  modes  designated  as  residual  show  very  little 
movement  as  a  result  of  any  of  the  control  forces  required  or 
transformations  applied  to  the  various  systems.  As  a  result 
of  the  choice  of  the  higher  frequency,  modes  as  residual  is 
verified . 


Tn t  roduc t i on 

With  the  success  of  the  Space  Shuttle  Program,  we  have 
entered  an  era  where  the  construction  of  large  space  struc¬ 
tures  will  become  a  reality.  To  achieve  practicality  and  use¬ 
ful  system  efficiencies,  the  proposed  sizes  of  these  struc¬ 
tures  are  hundreds  of  meters  in  diameter.  As  the  size  and 
flexibility  of  these  structures  increase,  the  number  of  low 
frequency  structural  modes  that  enter  the  bandwidth  of  system 
controllers  also  increases.  To  accomplish  control  of  such 
vehicles,  modeling  becomes  very  critical.  Even  with  improved 
modeling  techniques,  there  are  still  modeling  inaccuracies 
which,  in  the  limit,  could  result  in  unstable  conditions  if 
not  properly  compensated. 

The  method  of  control  that  is  both  realizable  and 
viable  is  modern  state  space  control  theory.  Using  this 
method,  however,  due  to  computat ; onal  requirements,  only  a 
Limited  number  of  structural  modes  can  be  handled  by  any 
single  controller.  Hence,  reduced  order  controllers  are 
required.  The  coupling  of  these  reduced  order  controllers 
with  detailed  finite  element  analysis  of  the  particular 
structure  can  be  successfully  adapted  to  meet  the  require¬ 
ments  of  several  missions  and  varied  flexible  structures. 

The  limiting  factor,  as  to  how  many  of  the  finite 


number  of  modeled  modes  may  be  successfully  controlled,  is 


the  capabilities  of  the  on-buard  tuii;|)iii  i-r  .  As  i  r>  ■  j  ]  i  •>: 
these  limitations,  only  these  modes  which  are  deemed  detri¬ 
mental  to  mission  requirements  are  controlled.  A  specific 
example  would  he  a  photographic  satellite-  where  point. ins 
accuracy  is  considered  critical  while  minor  vertical  vibra¬ 
tions  may  be  considered  inconsequential;  as  a  result:,  only 
those  modes  affecting  pointing  accuracy  would  be  controlled. 

While  specific  control  of  these  isolated  modes  would 
be  ideal,  it  must  be  realized  that  in  the  real  situation,  the 
sensor  data  will  be  contaminated  by  the  uncontrolled  modes 
and  these  same  uncontrolled  modes  may  be  affected  by  required 
inputs  to  the  desired  modes.  These  coupling  affects  are  refer¬ 
red  to  by  Balas  (Ref  1)  as  "observation  spillover"  and  "con¬ 
trol  spillover".  It  is  shown  that  either  of  these  system 
coupling  effects  may  lead  to  overall  system  instability.  The 
method  of  control  proposed  by  Balas  is  based  on  the  use  of 
narrow  bandpass  filters  which  effectively  comb  out  the 
suppressed  modes,  thus  eliminating  observation  spillover. 

Another  technique  which  was  first  presented  by  Sesak 
(Ref  2)  and  later  expanded  on  by  Coradetti  (Ref  3)  involves 
a  "singular  perturbation"  technique.  It  is  concluded  that 
this  approach,  with  infinite  penalty  on  spillover,  is  essen¬ 
tially  the  same  as  finding  transformation  matrix.  By  apply¬ 


ing  this  transformation  matrix  to  the  associated  gain 
matricies,  either  controller  or  observer,  the  spillover  terms 


would  bo  driven  to  zero.  This  method  can  be  effective  in 
removin',  destabilizing  cross  coupling  terms  even  if  these 
terms  do  not  result  in  overall  system  instability;  thus 
improving  system  response.  These  goals  are  accomplished 
through  the  application  of  state  space  control  techniques 
coupled  with  singular  value  decomposition  of  rectangular 
matrices  of  modal  amplitude  (Ref  4). 

The  primary  thrust  of  this  investigation  is  to  study 
the  application  of  the  above  techniques  on  the  implementation 
of  two  or  more  decentralized  controllers  on  a  lumped  mass 
model  of  a  tetrahedron.  The  primary  means  of  evaluating  the 
effectiveness  of  the  system  will  be  an  eigenvalue  analysis 
of  the  closed  loop  system  and  a  time  response  of  the  pointing 
angles  to  initial  conditions.  This  work  first  investigates 
all  of  the  results  of  Janiszewski  (Ref  5)  and  then  expands 
from  the  single  controller  model  utilizing  only  eight  modes 
to  the  multiple  controller  system  using  a  twelve  mode  model 
representation.  The  elimination  of  any  spillover  terms  will 
be  accomplished  through  the  implementation  of  the  transforma¬ 
tion  technicue  mentioned  earlier. 

The  specifics  of  the  system  model  used  in  this  investi 
gation  will  be  fully  explained  in  the  following  section.  The 
model  is  configured  with  sensor  actuator  pairs.  The  sensors 
are  position  sensors  only  and  are  used  to  determine  the  modal 
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amplitudes  at  a  point.  A  singular  value  decomposition  is 
performed  on  the  matrices  of  modal  amplitudes  to  obtain  a 
transformation  matrix  which  is  employed  to  eliminate  spill¬ 
over  terms.  With  the  addition  of  a  second  controller,  the 
improvement  of  response  in  the  structure  is  dramatic. 

Finally,  the  possibility  of  implementing  more  than  two  decen¬ 
tralized  controllers  is  examined. 
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System  Mode  1 

The  model  used  for  this  investigation  is  the  tetra¬ 
hedral  model  developed  at  the  Charles  Stark  Draper  Laboratory, 
Inc.  This  model  was  arrived  at  due  to  the  fact  that  it 
not  only  displayed  many  of  the  characteristic  responses 
observed  in  large  space  structures.  It  also  provided  a  low 
order  model  upon  which  various  control  systems  could  be  easily 
applied  so  evaluation  is  simple  as  a  result  of  the  small  num¬ 
ber  of  modes  present.  The  performance  criteria  of  the  model 
is  based  on  the  motion  of  the  structure  at  node  1.  This  is 
analogous  to  a  line  of  sight  evaluation  of  a  typical  optical 
system. 


The  finite  element  model  of  the  structure  is  displayed 
in  Figure  1.  The  structure  is  pin  connected  at  each  of  the 
nodes;  as  a  result,  it  is  only  capable  of  transmitting  axial 
forces.  A  Youngs  modulus  value  of  one  was  used  to  simplify 
the  stiffness  computation.  The  beams  are  considered  massless 
with  all  mass  located  at  nodes  1  through  4.  The  measured 
location  of  each  node  is  listed  in  Table  I. 

An  eigenvalue  anal>sis  of  tne  structure  was  accomplished 
on  the  NASTRAN  Computer  Program.  The  key  results  of  the 
analysis  are  listed  in  Table  II.  Tie  associated  eigenvetors 
are  listed  in  Appendix  A.  Table  II]  is  a  listing  of  the 
initial  conditions  that  were  applied  to  the  model  to  achieve 
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a  time  history  of  the  response.  For  the  purpose  of  this 
investigation,  it  is  assumed  that  i hese  values  are  applied  t< 
achieve  a  desired  pointing  requirement.  As  a  ’'esult  ,  all  of 
these  values  could  be  inputs  to  the  controller  prior  to  actu¬ 
ation,  thus  achieving  initial  conditions  on  all  of  the  error 
terms  of  zero.  The  development  of  these  error  terms  will  be 
covered  in  the  following  derivations  of  the  equations  of 
motion . 
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Table  I 

Node  Coordinates 


Mode 

X 

Y 

Z 

1 

0.0 

0.0 

10.165 

2 

-3.0 

-2.887 

2.0 

3 

5.0 

-2.887 

2.0 

4 

0 . 0 

5  .  77  35 

2.0 

5 

-6.0 

-1  .  1547 

0.0 

6 

-4.0 

-4 . 6188 

0 . 0 

7 

4.0 

-4.6188 

0.0 

8 

6.0 

-1.1547 

0.0 

9 

2.0 

5 . 7735 

0.0 

10 

-2 . 0 

5.7735 

0.0 

Table  II 

Kesu  1 1  .i 

of  NASTRAN  Analys 

i  s 

General i zed 

Generalized 

W  /RAD\ 

n (sec/ 

5. /RAD  \  2 

Mode 

Mass 

St i f  fness 

VSEC/ 

1 

1.0 

1 . 37 

1 .  171 

1 . 37 

2 

1.0 

2.15 

1.467 

2  .15 

3 

1.0 

8.79 

2.965 

8.79 

4 

1.0 

12.6 

3.558 

12.6 

5 

1.0 

14.8 

3.848 

14.8 

6 

1.0 

26.5 

5.149 

26.5 

7 

1.0 

32.2 

5.676 

32  .2 

8 

1.0 

32.6 

5.711 

32.6 

9 

1.0 

79.9 

8.940 

79 . 9 

10 

1.0 

106 

10.030 

106 

11 

1.0 

119 

10.923 

119 

12 

1.0 

195 

13.966 

195 
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Equations  of  Motion 

The  equations  of  motion  for  the  vibrational  motion 
of  a  large  space  structure  can  be  written  as: 

m£  +  E£  +  K£  =  Du  (1) 

where  g  is  an  n-vector  of  generalized  coordinates ,  M  is  an 
n  x  n  summetric  mass  matrix,  K  is  an  n  x  n  symmetric  stiffness 
matrix,  u  is  an  m-vector  of  inputs,  D  is  an  n  x  m  matrix  of 
modal  amplituder  evaluated  actuator  locations,  and  E  is  an 
n  x  n  dumping  matrix. 

Rewriting  equation  (1)  in  a  modal  coordinates 

?  +  2  ^  “  'l  +  J  7  =  0T  D  u  (2) 

where 

g.  =  0  q  (3) 

T 

and  0  is  the  transpose  of  the  n  x  n  model  matrix  for  equation 
(1) .  The  model  matrix  0  is  such  that 

g1  m  2.  =  [i] 

0T  K  0  =  [w?] 

0 T  E  0  =  [2C  wj 

where  all  matrices  which  are  displayed  are  n  x  n  diagonal. 

To  be  more  explicit,  flj  is  the  identity  matrix,  L^l  is  a  matrix 
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of  Che  eigenvalues  of  equation  (1)  andjj^c-jjis  the  associated 
damping  matrix. 


By  placing  equation  (2)  into  state  vector  format,  we 
arrive  at  equation  (4) : 


whe  re 


x  =  A  x  +  B  u 


(4) 


In  the  practical  case,  the  complete  state  vector  is 
not  available  and  quation  (4)  must  be  supplemented  by  an  out¬ 
put  equation.  If  we  assume  both  position  and  velcoity  infor¬ 
mation  is  available,  the  general  output  equation  becomes: 

Y  =  Cp  g  +  Cv  g  (5) 

in  state  vector  from 

Y  =  C  x  (6) 


where 


v 


0 


Equations  (4)  and  (6)  are  tie  model  of  the  satellite  avail¬ 
able  to  the  control  designer. 
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Reduced  Order  Model 


The  state  vector  x  from  above  is  a  2n-vector  that 


represents  the  entire  structural  model.  This  state  vector 
can  be  broken  into  a  number  of  more  specific  portions  in  the 


form 


-  T 

x  ,  x 
s  r 


-  T 
,  x 

urn  _ 


X  - 

c 

-*2 

n 

c 

-  vector, 

Control  led 

modes 

x  — 

s 

~f2 

n 

s 

-  vector , 

Suppressed 

modes 

X  r 

-*2 

n 

T 

-  vector. 

Residual  modes 

x  *—-*2  n  -  vector,  Unmodeled  modes 
urn  um 


The  unmodeled  modes  are  mentioned  here,  but  will  not 
appear  in  any  of  the  derivations.  These  are  modes  which  are 
beyond  the  capability  of  the  mathematical  system  model  to 
approximate.  The  unmodeled  modes,  while  they  may  exist  in 
nature,  have  not  been  taken  into  the  modeling  effort  of  the 
system.  The  response  of  the  actual  system  will  determine  if 
a  more  indepth  model  is  required  to  mathematically  approximate 
the  real  world  system.  The  controlled  modes  are  those  modes 
determined  to  require  active  control  to  achieve  desired  system 
response  characteristics.  In  general,  these  will  not  neces¬ 
sarily  be  the  lowest  frequency  modes. 


Due  to  computational  limitations  in  control  of  the 


vehicle,  only  a  small  subset  of  the  modeled  modes  may  be  con- 


trolled,  as  a  result  ,  spillover  due  to  the  remain  in;-  inodes 
will  occur.  To  eliminate  the  Hi  la tor  ions  effects  of  this 
spillover,  a  portion  of  the  modeled  but  uncontrolled  mode 
w. 11  be  suppressed.  The  number  of  which  can  be  suppresed  is 
again  dependent  on  hardware  limitations.  The  remaining 
modeled  modes  will  be  termed  residual  modes.  These  residual 
modes  will  have  "spillover"  terms  and  can  be  considered 
representative  of  those  higher  frequency  modes  that  are 
unmodeled . 

With  these  definitions,  the  system  state  may  be 
represented  by  the  following  equations 


X 

c 

=  A 

c 

X 

c 

+ 

B 

c 

U 

(8) 

X 

s 

=  A 

s 

X 

s 

+ 

B 

s 

u 

(9) 

xr 

=  Ar 

X 

r 

+ 

B 

r 

u 

(10) 

Y 

=  C 

c 

X 

c 

+ 

C  ^ 
s 

f  C  x 
s  r  r 

(ID 

The  system  parameters  matrices  A,  B,  and  C  are  as 
previously  described  in  the  text.  Furthermore,  the  matrices 

B  can  be  expressed  as:  B  =  --2- 

c  C  0  D 

c 


where  •. 


where  .  is  the  location  of  the  j  '  ^  actuator  and  n  is  the 
i  a 

number  of  actuators  and  nc  is  the  number  of  controlled  modes. 

In  this  representation,  the  coefficient  matrix,  D,  has  already 

been  multiplied  through.  Those  values  for  the  structure 

studied  are  the  NASTRAN  results  in  Appendix  A.  Similar  repre- 
T 

sentations  of  0  D  can  be  developed  for  both  the  suppressed 
and  residual  modes  with  the  only  difference  occurring  as  n^ 
and  n  ,  the  number  of  suppressed  and  residual  modes,  respec= 
tively  are  used  in  place  of  n  . 

In  a  similar  manner,  if  we  assume  only  point  sensors 
located  at  the  points  f, .  then 


where 


C 

c 


0 


Cp  0  "  |"al  (^1)  . anc  ('1} 

1  c 

al  ^'nsen^  . anc  ^nsen^ 

where  n  is  the  number  of  sensors  employed, 
sen  r 

The  equations  developed  to  this  point  are  quite 
general  and  independent  of  structural  complexity.  With 
increased  complexity,  the  sizes  of  the  respective  matrices 
are  the  only  variable;  t:  h  a  t  will  increase  in  dimension. 

The  general  nature  of  the  development  is  the  key  to 
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it;;  wide  area  of  possible  application  ■  n  a  variety  of  iarre 
space  structures. 

Modal  Control 


The  control  modal  upon  which  the  control  design  is  to 
be  based  is  given  by 


X  =  A  x  +  B  U 
c  c  c  c 


(8) 


Y  -  C  x  +  CX+CX 
c  c  s  s  r  r 


('ll) 


Three  feedback  controllers  are  examined.  The  form  of 
control  for  each  controller  is: 


A 

°1 

-  G 

X 

c 

(12) 

A 

A 

(13) 

u2 

=  G1 

X1 

+ 

G2  X2 

A 

A  A 

(14) 

g3 

=  G1 

X1 

4- 

G2  X2  +  G3  X3 

A 

A  A 

where  X^  , 

X2  ,  and 

a  re 

the 

specific  modes  controlled 

by  each 

controller.  Ideally  the  control  law  would  be  G  X  but  in  this 
case,  where  not  all  of  the  states  are  available,  the  estimated 
values  of  the  states  must  be  used.  The  individual  closed  loop 
system  matrices  will  be  developed  sequentially  in  the  follow¬ 
ing  discuss  ion . 


Single  Controller 


Since  we  are  unable  to  directly  measure  the  entire 


state  vector,  it  is  nect-ssar  7  to  et:i;,  ]  o  an  ohai  .a  r  of  the 


f  or  Ill  : 


whe  re 


A  X  +  B  r?  f  K  (■:  -  7) 

e  c  c 


A 

i;  x 

c  c 


e  s  t  i  ina  ted  s  t  a  t  e  ve  c  t  o  r 


estimated  output  vector 


The  observer  pain  matrix  is  chosen  such  that  the  error 
in  the  state  estimate,  represented  by 


e  =  X  -  X 
c  c  c 


is  asymptotically  stable 


The  closed  loop  system  stability,  including  controller 
and  observer,  can  be  evaluated  by  writing  the  state  equations 
for  an  augmented  state  vector  defined  below.  For  the  single 
controller  Z  will  be  defined  as  follows: 


T  -T  —  T  -  T 
(X  ,  e  ,  X  ,  X  ) 

c  s  r 


With  the  definition  the  overall  closet’  loop  system 
matrix  can  be  represented  in  block  matrix  form  as: 

k  +  Bcc;  BcG  0  ol  < 


A  -  KC  KC 

c  c  s 


0 


A.  this  point  it  is  of  interest  to  1  ool-  at  the  develop¬ 
ment  of  tlie  observer  sain  matrix,  K  and  the  control  feedback 
gain  matrix,  (’. .  First  consider  the  control,  sain  matrix  G. 


In  order  to  use  linear  optimal  regulator  theory,  a 
performance  index  is  defined  as: 


.1 

where 

Q 


1/2 


-  T 

X  o  X 
c  c 


+ 


R  u )  dr 


(20) 


is  an  n  x  n  positive  semide finite  weighting 
matrix. 


R  -  is  an  m  x  m  positive  definite  weighting 
matrix. 


This  performance  index,  subject  to 

X  =  A  X  +  B  L’ 
c  c  c  c 

is  minimized  with 

U  =  G  X 

c 


and 


- 1  T 

G  =  -R  B  S 
c 


(21) 


and  S  is  the  solution  to  t  le  matrix  Recatti  Equation. 


sa  +  ats-:;b  r-1btx+q=o 

C  C  C  C  x 


(22) 


The  development  of  the  observer  matrix  can  be  formu- 
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formulated  in  an  identical  development  once  it  is  realized 
that  the  eigenvalues  of  the  matrix,  (A  -  KC^) ,  are  equal  to 
the  eigenvalues  of  the  transpose  of  the  matrix.  The  system 
can  be  then  written  as: 

W  (t)  =  Acr  W  (t)  -  CT  g  (t) 

g  (t)  =  KT  W  (t) 

Using  this  system  and  defining  a  similar  performance 

index  as  listed  in  equation  (20)  with  the  substitution  of  W 

T 

for  Xc  leads  to  the  solution  for  the  gains  K  in  the  form. 

_ T  _  _  1  _  - 

K  =  +  R  ,  C  P 

ob  c 

where  P  is  the  solution  to  the  steady  state  algebraic  matrix 
Ricatti  Equation  .- 

pat  +  a  p-pctr,'1c  P  +  Q  ,  =  0 

c  c  c  ob  c  ob 

While  the  system  of  equations  is  not  block  triangular 
it  can  be  made  block  triangular  through  the  elimination  of 
control  spillover  or  observation  spillover.  Once  we  have 
achieved  suppression  of  the  appropriate  terms,  the  stability 
of  the  system  is  assured  through  the  proper  design  of  the  con 
troller  and  observer.  For  the  purpose  of  this  research, 
elimination  of  observation  spillover  has  been  deemed  more 
practical  and  cost  efficient.  Additional  sensors  to  achieve 
the  desired  observation  spillover  is  much  easier  to  implement 
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than  in  ere::  si  nit  the  number  of  ,tcf  u  1 1  u  in  achieve  spill¬ 
over  suppression. 

Dual  Controller 

The  following  development  of  a  two  controller  system 
parallels  that  of  the  single  controller.  The  control  law  to 
be  applied  is  as  stated  in  equation  (13).  In  this  system 
rather  than  defining  a  specific  number  of  modes  as  suppres¬ 
sed,  the  goal  is  to  achieve  two  decentralized  controllers 
which  will  be  independent  of  each  other. 

The  two  state  equations  are: 

XL  =  Ax  Xx  +  B1  U  (23) 

X2  =  A2  x7  +  B2  U  (24) 

Recalling  the  general  observer  equation  (13)  and  equa¬ 
tion  (16)  where  the  control  law  applied  is  equation  (13). 

X.  =  A.  X.  +  B.  U  +  K  (y  -  y)  i  =  1,  2 

l  it  t  J ' 

A  A 

Y .  =  C .  X .  i  =  1  ,  2 

tit 

_  A  A 

U  =  G1  Xx  +  G2  X2 

The  error  in  each  system  is  described  as 

-  X.  i  *  1,  2  (25) 
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By  apply  in 

a'  equat  ion 

a  (11) 

,  ( 

2  3) 

.  (2/.).  (2  3) 

,  and  t  he  i-s t  i ma tor 

equat .ons 

which  are 

1  i  s  t  e  d 

above 

,  i  -fc  can  be 

;hown  that  e.  is 

t 

described 

by  : 

/\ 

®1  = 

X1  -  Xi  = 

(A1  - 

K1 

C1 

) e i  +  k1c2x2 

+  K,C  X  (26) 

1  r  r  7 

*1  - 

x2  -  xL  = 

(A9  - 

k2 

C2 

+  K2C1X] 

f  K„C  X  (27) 

2  r  r 

The  associated  X  equation  may  be  simply  derived  using 
the  system  equation  (23)  and  the  control  law  (13).  The  result¬ 
ing  equation  is: 


Xx  =  (Aj_  +  Bj  Gl)  Xj_  +  Bx  G:  e1  +  B1  G?  c'2  +  B1  G?  X, 

(28) 


A  similar  application  of  the  control  law  and  the  residual  model 
equation  (10)  provides  the  following  results: 


X  = 
r 


A  Xr  +  B  G,  X,  +  B  G,  e,  +  B  G0  X„  + 
rr  rll  rll  r22 


B  C, 
r  Z 


(29) 


By  defining  an  overall  system  vector  z  of  the  fom: 


Z  r 


T  -  T  T  -T  T 

r 


X1  ’  el  -  X2  ’  e2  ’  X 


(30) 


The  closed  loop  system  model  including  the  two  decentralized 
controllers,  each  utilizing  state  variable  feedback,  can  be 
wr  i t ten  as  : 
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0 


(31) 


Z  = 


A1+  -h  Gl  !31  Gx 


Bl(:2  B1  C2 


B2  ('l 


K2  C1 


B  G, 
r  1 


(A1-K]C1)  K  L  C , 


b2gl  (a2+  b2g2)  b,  g2 


K.C 
1  r 


(A,.-K0Co)  K0C 

2  2  2  2  r 


B  G, 
r  1 


B  G„  B  G,, 

r  2  r  2 


Z(t) 


It  is  apparent  that  the  suppression  of  all  the  "obser¬ 
vation  spillover"  or  the  "control  spillover"  terms  is  insuf¬ 
ficient  to  completely  triangularize  the  system  even  in  the 
absence  of  residual  modes.  To  achieve  a  closed  loop  system 
with  the  above  characteristics,  it  is  necessary  to  suppress 
the  control  spillover  term  of  one  system,  e.g.,  B^  ,  while 
suppressing  the  observation  spillover  of  the  other  system, 

C2 -  The  judicious  selection  of  modes  again  is  critical  so 
as  to  provide  a  frequency  separation  between  the  lower  frequency 
controller  and  the  residual  modes. 


The  primary  advantage  of  two  controllers  is  the  number 
of  modes  controlled  can  be  divided  between  the  two  systems. 
This  is  important  since  the  computational  burden  of  solving 
the  Ricatti  Equation  increases  roughly  as  the  cube  of  the 
orcer  of  the  equation  (Ref  3).  Therefore,  the  advantages  of 
solving  the  Ricatti  Equation  of  two  smaller  controllers  is 
apparent . 


2 


Three  Con 1 1 olio r  s 


To  avoid  a  repetition  of  all  of  the  equations  developed 


i  n 

the  pre\ 

Lous  section 

,  it  can 

be  stated 

that  the 

control 

1  aw 

of 

equation 

(14)  was 

applied  to 

arrive  at 

the  closed  loop 

sys- 

tern  model  of  this  section.  The 

r  T  T  T  - 

Z  =  [Xx  ,  v1  .  X2  •  C2 

Z  ve dor 

T  ,,  T 
.  X,  .  e-5 

is  defined 

T-  Xr‘T 

(A  S  ! 

132) 

This  result 

s  in  an  ovei 

all  closed  loop  system  equation: 

At+BtC, 
l  i  1 

B1G1 

B1G2 

B1G2 

B 1 G  3 

B1G3 

0 

0 

A1~K1C1 

K1C2 

0 

Kl°3 

0 

K,C 

1  r 

b2Gi 

B2G1 

A2+B2G2 

B2G2 

B2G3 

B2G3 

0 

K2C1 

0 

0 

A  2 -K9C2 

K2C3 

0 

K0C 

2  r 

B3G1 

B3G1 

B3G2 

B3G2 

A3B3G3 

B3G3 

0 

K3C1 

0 

0 

C 

0 

A3-K3C3 

K3cr 

BrG  1 

B  G, 
r  1 

B  G~ 

r  2 

B  G0 
r  3 

B  G. 
r  3 

BrG3 

A 

r 

As 

discussed 

earlier,  the  system 

presented 

here  cannot 

be  trianp;ularized  ti.rou  ’h  complete  elimination  ol  observation 
spillover  or  control  spillover.  There  are  two  approaches  that 
can  be  utilized  in  the  xamination  of  the  three  controller 


system.  First,  through  the  judicious  positioning,  of  sensors, 
the  modal  amplitude  matrix  and  thus  the  system  parameter 


(33) 


2  2 


matrices,  ft  and  C,  of  on  controller  can  be  made  orthogonal  to 


the  remaining  two  controllers.  To  completely  decouple  the  sys 
tern,  the  terms  which  must  be  eliminated  are  listed  in  Table 
IV.  By  arranging  the  modes  such  that  two  of  the  controllers 
operate  on  modes  such  that  two  of  the  controllers  operate  on 
modes  that  are  orthogonal  to  the  third  controller,  the  system 
would  reduce  to  a  two  controller  system.  The  system  model 
used  in  this  study  has  been  determined  to  contain  such  proper¬ 
ties.  This  will  be  specifically  demonstrated  in  the  investi¬ 
gation  portion  of  the  text.  At  this  point  suffice  it  to  say 
that  the  twelve  modes  modeled  can  be  divided  into  two  ortho¬ 
gonal  groupings. 

As  an  example,  let  controllers  one  and  two  operate  on 
the  first  group  of  orthogonal  modes  while  controller  two 
operates  on  a  portion  of  the  second  grouping.  As  a  result, 
all  cross  terms  between  one  or  two  and  three  will  be  equal  to 
zero.  This  will  reduce  Table  IV  to 


B1G2 

=  0 

E2G1 

=  0 

or 

K1C2 

-  0 

K2C 1 

=  0 

whicl  are  the  terms  required  equal  to  zero  to  decouple  the 
two  controller  system,  therefore  demonstrating  the  ability 
to  reduce  the  system  to  a  two  controller  problem.  The  second 
method  of  system  suppression  would  require  an  optimization 
process  included  in  the  tr msformat  i  01  formation  such  that 
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- 


Table  IV 


Total 

Decouple  of  3  Controller 

B2G3 

=  0 

B2G1 

K2C3 

=  0 

k2c] 

B1G3 

=  0 

B3G1 

OR 

K1C3 

=  0 

k3c1 

B1G2 

=  0 

B3G2 

K1C2 

=  0 

k3c2 
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such  terms  us  B7Gj  and  B  Uj  <lrr  approximately  zero.  While  the 
possibility  of  obtainin'’  a  transformation  matrix  orthogonal  to 
both  matrices  is  highly  unlikely,  an  opt i mazat i on  process  can 
be  applied  to  reduce  the  value  of  these  spillover  terms  to 
insignificant  values  relative  to  the  system  dynamics.  This 
second  nethod  would  require  more  on-board  computational  capa¬ 
bilities  which  may  result  in  exceeding  the  designed  capacity 
of  the  system.  As  a  result,  this  method  would  be  far  more 
costly  to  implement  thus  making  the  first  method  the  only 
viable  approach.  Since  it  has  been  demonstrated  that  the 
system  can  always  be  reduced  to  a  two  controller  problem 
through  proper  sensor  and  actuator  location,  only  the  investi¬ 
gation  of  the  single  and  dual  controllers  will  be  carried  out 
in  this  research. 
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1 


Trans  format i on  Technique 


This  section  is  designed  to  describe  in  further 
detail  those  methods  applied  to  the  model  to  achieve  the 
block  triangular  form.  This  will  require  the  elimination  of 
the  cross  coupling  terms  such  as  and  B  j  G  ^  in  the  two 

controller  system.  The  entire  thrust  of  this  method  is  to 
drive  these  terms  to  zero  while  keeping  the  terms  B-^G^  , 

B2G2  ,  and  ,  K.2C2  not  equal  to  zero.  This  is  first  done 

for  the  single  controller  case.  The  technique  is  then  applied 
to  the  two  controller  problem  by  eliminating  the  control 
spillover  of  one  controller  while  operating  on  the  observa¬ 
tion  spillover  of  the  second  controller. 


For  the  single  controller  system,  the  elimination  of 
observation  spillover  is  achieved  if  a  R  matrix  can  be  found 
such  that 

K  C  =0  (34) 

s 

K  C  =0  (35) 

r 


while 


K  C  =  0 

c 


(36) 


The  final  equation  is  constraint  that  must  be  met  in 
order  to  maintain  observability  over  the  controlled  modes. 

While  it  would  be  optimal  to  achieve  both  equation 
(34)  and  equation  (35)  in  the  system  model,  in  the  actual 
structure  thif  would  not  be  fully  realizable.  This  is  primarily 
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due  Co  the  large  number  of  modes  that  are  ;.h  v,  i  iu !  !  .•  present 
in  the  structural  model.  As  a  result  ,  on  1 7  a  subset  .  >  i  tin- 
modeled  modes  will  be  suppressed.  Thus,  on!-/  equal  i  on  {  , 

will  be  sat  is f ied . 

The  selection  of  those  inode.1  to  be  designated  as 
suppressed  or  as  residual  is  somewh;  t  arbitrary  and  could  be 
established  through  an  iterative  pr<  cess  .  Those  modes  you 
are  most  interested  in  suppressing  ere  those  modes  which, 
even  though  stable,  are  weakly  damped  and  thus  may  be  driven 
unstable  as  a  result  of  the  observation  spillover.  The 
selection  of  those  modes  as  residual  would  be  best  designated 
as  those  modes  which  are  actually  shifted  further  to  the  left 
of  the  jw-axis  as  a  result  of  the  observation  spillover,  thus 
stabilizing  these  modes.  Another  choice,  the  one  used  in  this 
investigation,  is  to  suppress  all  uncontrolled  modes  below  a 
certain  frequency.  The  primary  assumption  here  is  that  the 
higher  frequency  modes  fall  outside  of  the  bandwidth  of  the 
controller . 

However  the  selection  of  the  suppressed  modes  is 
accomplished,  the  system  to  be  examined  is: 

C  T  KT  -  0  (37) 

) 

This  is  nothing  more  than  the  transpose  of  equation  (34), 
however,  this  fori:  f  equation  is  more  useful  as  will  become 


apparen t 


To  achieve  this  do  si  red  result.,  the  K 


(■'jua  - 


ma  I  r  i o  t 

tion  (  7)  must  be  transformed  such  that  it  is  orthogonal  to 
T  T 

the  .\y  s  of  C  (columns  of  C  )  .  The  C  matrix  is  sized  such 
s  s  s 

that  it  has  the  number  of  columns  that  corresponds  to  the  num¬ 
ber  if  sensors  (n  )  and  a  non  zero  row  length  of  the  number 
of  s  ip ;  ressed  modes  (n  )  .  Looking  at  the  equation  of  the 
transformation  required 

CsTt  =  0  (38) 

The  number  of  linearly  independent  algebraic  solutions,  t, 

T 

are  :,p<  cified  as  the  difference  between  the  rank  of  and 

nsen  The  number  of  suppressed  modes  is  equal  to  or  greater 
than  the  number  of  sensors,  no  solution  vector  t  can  be  found 
unless  the  rows  of  are  not  linearly  independent.  As  a 

result  of  this  relation,  in  general,  the  number  of  modes  that 
can  :ie  suppressed  can  not  exceed  the  number  of  sensors  avail¬ 
able.  In  terms  of  output  we  will  define  v  by: 


As  a  result  of  the  out  out  •.  does  no  (  contain  the  sup¬ 
pressed  modes.  The  new  control  problem  to  be  eons  i  tic  ad  can 
be  stated  as : 

X  -  A  X  f  B  C  (U>) 

c  c  c  e 

v  =  T  C  X  +  T  C  X  I--  X  +  C"'  X  (42b) 

c  c  r  r  c  e  r  r 

The  output  v  is  no  longer  a  /ector  of  dimension  n 

T 

but  has  dimension  (nsen~  Rank  of  Ct,  >.  The  suppression  may, 
therefore,  be  thought  of  as  replacin',  a  system  of  n  sensors 
with  nsen'r)  synthetic  sensors. 

As  long  as  the  system  of  equations  (42)  are  observable 

and  controllable,  the  stable  matrice  .  A  +  B  G  and  A  -K  C 

c  c  c  c 

can  be  formed  and  placed  in  the  over  ill  system  matrix  of  equa¬ 
tion  (19)  in  which  the  observation  snillover  will  have  been 
removed.  If  the  suppressed  modes  for  this  system  are  properly 
chosen,  the  entire  system  will  remain  stable. 

With  that  general  overview  of  the  purpose  an  re '•nit 
of  the  technique,  the  specifics  of  obtaining  the  matrix  T  will 
be  developed.  The  matrix  of  interest  in  this  technique  v.  C 
observation  suppression  is  C  .  This  matrix  can  be  written  in 
the  form: 

C  =  W  t  vT 
s 

2  9 


(43) 


where : 


and 


W  i  S  an  n  ri 

.  i  n  in 

urt  hoj'.on.'i  1  matrix  of  le i  t 
s  i  mm  1  ar  vec  t  or:;  . 

V  is  an  n^  x  n_  orthogonal 
matrix  of  r i gh t  singular 
vectors . 


(V+) 


such  that 

of  C  and 

s 

more  ,  the 


S  m  a  r  x  r  matrix  of  the  non-zero  singular  values 
r  i  ;  t lie  rank  of  Ct.  as  previous!;/  stated.  i'urther- 
matrix  W  can  be  partitioned  such  that: 


W 


W  ;  W 


L r 

The  partition  W  is  an  n  x  r  matrix  of  left  si 
r  sen 

lar  vectors  associated  with  the  non-xero  singular  values 
and  Wq  is  an  nsen  x  q  matrix  of  left  singular  vectors 
associated  with  the  zero  singular  values. 


(45) 

ngu- 

o  f 


and  W 

r 


, .  T 

Since  W  is  an  orthogonal  matrix,  the  product  of  W i 

Ci 

is  zero  which  leads  to 


i  s 


W 


T 


C 


q  S 

As  a  result,  it 
composed  of  the  left 


T  T 

=  W  W  S  V  =  0 

q  r  r 

is  obvious  that  the  T  matrix  sought 
singular  vectors  associated  with  the 
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Computer  Model 


Th.'  primary  goals  in  the  formulation  of  the  program 
wer  ;  flexi  liLity  and  simplicity.  The  program  in  capable  of 
making  several  diverse  runs  depending  on  Lhe  desired  output 
or  the  particular  area  of  interest  being  examined.  The  pro¬ 
gram  generates  output  data  for  a  single  controller  or  a  dual 
controller  model.  In  either  of  these  types  of  runs,  the 
inclusion  of  the  residual  modes  is  optional. 

The  repetitive  nature  of  the  formation  of  the  para¬ 
meter  matrices  (A,  B,  C,);  in  the  control,  suppressed,  and 
residual  form;  led  to  these  being  structured  in  a  subroutine. 
This  format  also  added  the  flexibility  to  change  the  sixe  of 
the  matrices  to  meet  specific  requirements  of  various  investi¬ 
gations.  The  formation  of  the  initial  condition  vector  is 
also  accomplished  in  a  subroutine.  There  is  a  separate  sub¬ 
routine  for  the  different  initial  condition  vectors  required 
for  the  single  or  dual  controller  model. 

For  program  initialization  certain  data  is  required 
from  either  permanent  files  or  parameter  assignments.  Required 
data  from  the  user  is  the  number  of  controlled,  suppressed, 
and  residual  modes  followed  by  the  number  of  actuators  and 
sensors.  Finally,  the  damping  ratio  for  each  of  the  modes, 
assumed  equal,  must  be  designated;  0.005  for  the  system 
studied.  The  system  will  then  read  from  permanent  file  the 
NAS  'RAN  values  of  the  modal  amplitude  at  each  actuator  location 

r? 


Tilt-  same  data  is  then  entered  fur  each  of  the  sensor 
locations.  For  the  model  considered  these  values  are  identi¬ 
cal  since  the  sensors  and  actuators  are  colocated,  however, 
it  was  considered  necessary  to  make  separate  entires  to 
accommodate  those  possible  situations  where  the  sensors  and 
actuators  are  not  colocated.  Finally,  the  associated  frequency 
for  each  of  the  modes  is  read  in  from  permanent  file.  For 
time  response  calculation  the  initial  conditions  for  the  sys¬ 
tem  and  the  mode  shapes  of  the  point  of  interest  msut  also  be 
made  available. 

With  the  preload  of  this  data,  the  modal  arrangement 
as  controlled,  suppressed  or  residual  is  at  the  option  of  the 
operator.  The  various  modes  may  be  moved  in  any  manner  desired 
by  the  operator  without  the  requirement  for  the  preload  of 
any  additional  modal  information.  Once  a  particular  selection 
is  made,  the  program  will  form  the  specified  matrices  and 
the  associated  initial  condition  vector.  With  the  system  now 
completely  structured,  the  steady  state  feedback  matrices  are 
formed  (G  and  K) .  This  is  accomplished  through  the  execution 
of  a  series  of  sophisticated  subroutines  created  by  Kleinman 
Ref  6),  which  provide  a  numerical  solution  to  the  matrix 
Ricatti  Equation.  With  these  matrices  formed,  an  overall 
system  matrix  as  depicted  in  either  equation  (19)  or  equation 
(31)  is  formed  dependent  on  whether  a  single  or  dual  control¬ 
ler  run  has  been  indicated. 
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At  this  point  one  can  execute  the  option  to  create 
a  time  history  of  the  line  of  sight  a'  point  1  in  the  /.  and 
y  directions.  This  type  of  system  response  was  used  .since 


the  pointing  accuracy  of  the  vehicle  was  a  criteria  for  deter¬ 
mining  the  success  of  the  system  controller(s) .  The  line  of 
sight  was  calculated  with  the  use  of  the  zero  input  equation 
for  the  state  equation; 


A  x  +  B  u 

(46) 

C  X 

(47) 

The  zero  input  equation  is 


X  (kdt)  =  eAdt  X  (k-l)dt 


08) 


To  minimize  any  problems  that  might  arise  as  a  result 

of  rapid  system  osciallations  not  perceived  by  the  discrete 

A  d  t 

model,  a  DT  -  .01  was  utilized.  The  e  matrix  is  determined 
through  the  Taylor  series  expansion  of  the  term  in  the  fol¬ 
lowing  equation: 


Adt 


2  2 

I  +  Adt  +  A  dt  +  AJdt  + 


2  ! 


3! 


(49) 


The  value  of  the  displacement  at  position  1,  in  the 
x  and  y  direction,  is  calculated  each  .5  second  up  to  20 
seconds  using  the  m jde  shapes,  previously  loaded,  and  the 
computed  value  of  X  (t).  This  displacement  is  calculated 
through  the  summation  formula: 


Ill 


1  .  2 


ao) 


X  (t) 
n 


r 

irf 


n  1 


;  (t  ) 


Where  m  is  the  number  of  modes  in  model.  For  n  -  1 ,  the  equa¬ 

tion  computes  the  line  of  sight  displacement,  in  the  X  direction 
and  n  =  2  represents  the  Y  direction  displacement. 


At  this  point,  the  eigenvalues  of  the  matrices  (AFBG)  , 
(A-KC)  and  the  overall  system  matrix  are  computed.  This 
analysis  was  accomplished  by  implementing,  the  subroitine  EIGRF 
from  the  International  Mathematical  and  Statistical  Library 
(IMSL).  The  eigenvalues  of  the  closed  loop  system  matrix, 
since  they  reflect  the  overall  systems  stability,  determining 
success  of  the  system  suppression.  Comparison  of  t.iese  eigen¬ 
values  with  the  eigenvalues  of  (A+BC)  and  (A-KC)  demonstrate 
which  modes  were  most  affected. 

The  suppression  of  the  system  varies  whether  one  or 

two  controllers  are  implemented  by  the  operator.  With  a 

single  controller,  the  observation  spillover  is  eliminated  by 

T 

accomplishing  a  singular  value  decomposition  of  the  Cg 
matrix.  This  is  achieved  through  the  execution  of  the  IMSL 
subroutine  LSVDF .  Then  by  using  the  associated  singular  vec¬ 
tors,  a  transformation  matrix  is  generated.  The  overall  sys¬ 
tem  is  then  recreated  using  the  transformation  technique  as 
is  described  in  section  IV.  Once  the  new  gain  matrix  is 
created,  the  program  loops  back  and  initiates  another  time 


15 


response  and  eigenvalue  analysis. 

In  the  ease  of  the  two  controller  system  of  equation 
(31),  it  is  apparent  that  the  elimination  of  observation 
spillover  is  insufficient  to  completely  decouple  the  system 
in  the  absence  of  the  residual  mode: .  In  this  case  one  must 
eliminate  control  spillover  of  one  system  while  eliminating 
the  observation  spillover  of  the  other  system.  The  elimina¬ 
tion  of  system  one  observation  spillover  was  implemented  to 
take  advantage  of  that  part  of  the  computer  model  already 
structured.  The  control  spillover  are  implemented  in  the 
overall  system  and  a  time  response  calculation  followed  by  an 
eigenvalue  analysis  is  then  accomplished. 

Since  there  exists  the  possib  lity  that  modes  other 
than  suppressed  modes  are  adversely  affected,  a  series  of 
calculations  are  required  to  insure  tlat  controlled  modes  are 
not  aligned  with  any  of  suppressed  nodes.  A  further  check  to 
insure  that  controlled  modes  are  a  I  inear  combination  of  the 
suppressed  modes  is  accomplished  at  the  end  of  the  run.  When 
one  of  these  cases  are  encountered,  a  regrouping  of  the  modes 
is  required  to  avoid  the  detrimental  affects  of  the  suppres¬ 
sion  on  the  overall  system  response. 

While  the  model  used  in  this  study  is  very  specific 
in  its  definition,  the  subroutine  structure  of  the  program 
provides  the  flexibility  to  analyze  a  variety  of  other  sys¬ 
tem  through  the  simple  restructuring,  of  the  subroutines  that 


36 


t  e rn  to 


form  th<  A.  H ,  and  C  matrices  to  compl  /  with  t  hi-  new  sy:, 
be  anal-,  zed.  The  only  requiremei  t:  is  that  the  system  can  be 
written  in  the  format: 

X  =  A  x  +  B  m 

ant: 

,  Y  =  C  x 


Inves  tij^a  t.  i  on 

A  building  block  method  of  research  y;n>  deemed  t  he 
best  approach  to  make  a  thorough  study  of  the  complete  system. 
As  earlier  explained,  it  was  determined  that  due  to  hardware 
and  cost  considerations,  observation  spillover  elimination 
would  be  employed  when  at  all  possible.  Initially,  the  basic- 
system  was  researched  using,  only  controlled  and  suppressed 
modes.  This  was  done  to  confirm  the  fact  t  hat.  the  system  could 
be  block  triangularix.-d  through  the  elimination  of  observa¬ 
tion  spillover.  The  eigenvalue  analysis  that  resulted  is 
displayed  in  Table  V.  This  analysis,  wa •  accomplished  on  the 
first  eight  modes  of  the  NASTRAA st ructura 1  analysis.  The  con¬ 
trol  weighting  matrix  was  equal  to  Q  >()  |ij  . 

At  this  point,  the  four  residual  modes  were  added  to 
complete  the  imp  lemon  ation  of  the  twelve  mode  model.  The 
system  wa;;  analyzed  t  ;  determine  the  ffect  of  the  four  addi¬ 
tional  modes  on  the  response  characte- -i  st  i  cs  of  the  system. 

The  eigenvalue  analysis  of  this  syste  .  i listed  in  Table  VI 
with  the  associated  time  response  Ii  scd  in  Table  Vi  » . 

The  value  of  0  wax  left  at  tv-  value  of  .10  ^1  j  for 

the  remainder  of  the  investigation  :,<•  that  all  results  could 
be  associated  with  either  the  number  >f  controllers  imple¬ 
mented  or  the  variou  groupings  of  tie  nodes 
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Tab  1  c 


■  vs 

tf  .! 

12 i  ;;u!iv.il  tie 

Analysis  -  Single  Contra 

1  ;e 

r 

Moil  a  1  A:,  s  i  /nmen  i 

C( 

n  t ru 1  1 

eJ 

Suppressed 

s i dua 1 

1 ,2,4, 

5 

3 . 6 , 7 . 8 

lone 

Ove  ra  1 1 

System  Eigenvalues 

I'k 

■  ore  T 

rans 

f  o  rma  L  i  on 

After  Tran.s 

f  orma t i on 

- 

i007n 

f 

5 . 73519 

-  .028  3  7 

-t 

5 . 76583 

>  1 99  b 

+ 

5  .  75026 

-  .02855 

9 

5.71073 

12574 

9 

5 . 14935 

-  .  02574 

9 

5 . 14935 
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9 
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-  1  .  04  342 

9 

3 . 42286 

1 . 

12459 

+ 

3  .  27924 

-  .  34269 

9 

1 . 16273 

5796 

9 

1.25179 

-  .01482 

9 

2.96457 

■  5138 

+ 

0 . 77807 

-  .  50461 

9 

1 . 43353 

1 . 

'3795 

9 

3.76850 

-  .  13144 

9 

1.46583 

1 . 

38709 

9 

3 . 54688 
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9 

1 . 16990 
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9 

3.00503 

-1.02693 

9 
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+ 

1 . 58838 

-1  .24619 

9 

3.66533 
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9 

1.00300 

-1 .26459 

9 
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Table  Va 


Time  Response  -  Sinp.le  Controller 


Cont  ro 1  led 
1 .2  ,4,5 


Model  1  Ass  i  pnnien  t  s 

Suppressed 

3 .6 ,7 ,8 


Res i dua 1 
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Re  fore  Transformation  After  Trans ‘"ormat  ion 
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.000080 
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. 000408 

3.0 
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. 000641 
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.000533 

.000541 

4  •  5 

/  .  ,  1 

.  V  S  /J  i 

.  000633 

5-0 

.000327 

.000243 

p.  0 

.  ^08056 

. 00066  6 

5-5 

- .000408 

- .000874 

V  ■  ■' 

'  '•91,  V 

-  .  0002  57 

6.0 

- . 000202 

- .000622 

0 . 0 

-  •  O  _.-r- 

-  . 000269 

6.5 

. 000S44 

- .000322 

6 . 5 

-  .0  001-1 
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8.0 

-.000425 

.  000081 

8 . 0 

-  .  0  _ 

- .  000309 

8.5 
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.000614 

6 . 5 

.  000201 

9.0 

.000367 

.000485 

9-0 

J  'j  J  ' 

.  0002  92 

9.5 

- .000394 

.000082 

9  •  5 

- . 0000 3 l 

. 6 G 02 01 

10.0 

- .000469 

- . 000474 

10.0 

-  »  Owo  '/i  ) 

-  .  000354 

10.5 

.000206 

.000047 

10.5 

- . 00006 ‘ 

.  000084 

11.0 

.00063 

-.000112 

11.0 

.000^  57 

.  0 0 0 172 

11.5 

.000009 

- . 0001 09 

11.5 

.  JO 01  /  / 

.  00012s- 

12.0 

- .000402 

- . 000474 

12.0 

—  •  y 

-.000335 

12.5 

.000077 

.0001 57 

12.5 

-  .  0001  ?'■ 

-  .  Go 01 57 

13.0 

.000379 

.OCO329 

13.0 

.  uOo  j  ^ 

.0002 52 

13.5 

- .000034 

. 000182 

13.3 

r  f  *  1 

•  o/uO  y  + 

. 000123 

14.0 

- .000583 

- .0001 1 9 

14.0 
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14.  5 

- .000083 

- . 000046 

14.5 

-  •  0O0 

- . 000287 

15.0 

.000376 

.000317 

15.0 

.  0002  r*- 

. 0C0270 

15-5 

. 000226 

-  .0001 1 5 

1  r  ‘ 

.000393 

.  -  j 001 42 

10.0 

- .OOO358 

-  .000186 

16.0 

- .000186 

-  .  0  0  0  0  3  2 

16.5 

-.OOOO96 

-.000332 

16.5 

-  .  000463 

-  . OOO320 

17.0 

.000362 

.0002 99 

17-0 

.000167 

.000111 

17.5 

.000217 

.000008 

17.5 

. 00040  5 

.  0002  36 

18.0 

- .000953 

-  .  000026 

16.0 

. 000007 

-  . 00038 

1  ■ .  5 

- . 000322 

-  .000156 

18.5 

- . 000494 

-  .  00022  6 

I9.O 

.000263 

.  000251 

19.0 

. 000004 

-  .  000076 

19.5 

.000224 

. 000195 

1 9 . 5 

.000405 

.000326 

2e  .0 

-.000162 

-  .000166 

20.0 

.  00  01 2  9 

-  .  000014 
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Bast'd  on  the  aft  net  of  the  'undos  on  tin.-  motion  of  the 
structure,  it  was  deemed  most  beneficial  to  control  modes  1 


2  ,  4  ,  5 

wh  i  1  e 

suppres sing 

modes 

1  , 

6  ,  7 

,  8.  1 

■'  i  n  a  1  1  y  ,  the 

res i dun 

1  modes 

were  9  ,  10 

,  11  , 

12  . 

The 

eho i  ci 

o f  the  re s 

modes  was  based  on  the  fact  that  because  of  freqency  separa¬ 
tion,  these  modes  would  be  unaffected  by  the  control  that  was 
applied  to  lower  frequency  modes.  This  premise  is  born  out 
in  the  analysis  of  the  eigenvalues  presented  in  Table  VI. 

This  shows  a  damping  ratio  for  the  residual  modes  of  approxi¬ 
mately  .001  or  greater,  since  0.005  was  used  as  the  open  loop 
dampinp  ratio,  the  controller  has  increased  the  damping  of 
each  of  the  residual  modes  even  through  they  were  not  included 
in  the  optimal  control  formulation. 
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Taf  1 e  VI 
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Table  Via 


Time  Response 
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- .00001 7 

■o 


Values  are  shown  before  and  ifler  suppression 

The  next  logical  step  was  to  exami  no  the  system  per¬ 
formance  with  the  dual  controller  svstem  of  equation  ('ll) 
was  implemented.  Again,  this  research  was  done  with  0  -  2') 
Tl]-  Tiie  system  was  divided  such  that  controller  one  handled 
modes  1,  2,  4,  and  5,  as  determined  necessary  to  achieve 
acceptable  pointing  accuracies,  controller  two  was  initially 
specified  as  modes  3,  6,  7,  and  8.  The  two  controller  model 
was  run  agains  the  eight  mode  truncated  model  to  confirm  the 
effectiveness  of  the  method  of  suppression  employed.  These 
results  can  he  seen  in  Table  VII. 

With  these  results  the  aiditional  residual  mode  (9, 

10,  11,  12)  were  included  in  the  model  to  check  for  any 

adverse  effects  due  to  these  modes  as  was  encountered  in  the 
earlier  investigation  of  the  single  controller.  In  this 
case,  the  overall  system  retained  the  achieved  stability  as 
is  seen  in  Table  VIII.  The  desirable  results  that  were 
achieved  in  this  arrangement  were  that  the  system  achieved 
the  desired  accuracies  in  the  x  and  y  directions  within 
approximately  10  seconds.  The  associated  time  response 
printout  to  each  of  the  above  runs  are  presented  in  Tables 
IV^through  VIII*- 
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Appendix  A.  Using  the  definition  of  the  dot 
vectors  : 
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The  angles  between  t ne  modal  amplitude  vectors  were  deter¬ 
mined.  This  was  done  to  determine  if  any  of  the  modes  lied 
on  lines  of  action  such  that  they  could  either  be  simply 
separated  or  arranged  to  minimize  control  efforts  required 
by  associating  similarly  aligned  modes.  As  a  result  of 
this  inves t i ga t ion ,  it  became  evident  that  the  modal  ampli¬ 
tude  vectors  subdivided  into  two  orthogonal  vectors  ("Table 
9)  . 


With  these  orthogonal  groupings,  the  system  was  run 
with  controller  one  operating  on  modes  1,  A,  6,  and  7, 
while  controller  number  two  drove  modes  2,  3,  5,  and  8.  This 
grouping  provided  the  best  overall  system  response.  The 
eigenvalues  of  this  system  is  depicted  in  Table  X  while  the 
associated  time  response  is  listed  in  Table  Xa . 

The  unique  quality  of  this  system  is  that  it  is  in¬ 
herently  decoupled,  in  that  the  associated  feedback  gain 
matrices  of  one  system  (K  and  G)  are  orthogonal  to  the  other 
system  parameter  matrices  (B  and  C) .  This  results  in  the 
fact  the  off  diagonal  coupling  terms  G^  ,  K?  C2  , 

G2  are  all  equal  m  zero. 
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This  nodal  vector  ort  hogona  1  i  t.  •/  while 
occurancc ,  shows  the  impot [ ance  of  proper  location  of  the 
sensors  and  actuators  on  the  model .  This  system  analysis  to 
locate  the  sensors  on  the  structure  is  a  design  tool  which 
should  not  be  taken  lightly.  The  judicious  location  of 
sensors  and  actuators  can  reduce  the  system  to  a  pair  of 
uncoupled  controllers  requiring  no  system  suppression;  as  a 
result,  no  degradation  in  the  system  response  from  the  optimal 
gain  values.  By  referring  again  to  Table  7a,  it  is  obvious 
that  the  time  response  of  the  system  before  suppression  is 
superior  to  that  after  suppression. 


50 


,  4> 


S  /.■ 

.  t  e 

m  li  i  yen  v  t  1  m  •  An  i  ' 

.is  -  Twi  >  Control 

!  e  r.s 

*1  o  1 1  a  1  As:,  i 

■  nt'ien  1  s 

Cont  ro  l  Li 

■  r 

il  1  Control] 

or  2 

Ret; 

i  d  u  a  1 

1,4.6. 

7 

2  .  3  , ' 

>  H 

9.10 

,11,12 

Ove  r;i  11  Sy  s  L  cm 

0  i : ’,en  va  1  mcs 

Before  Tran 

s  format  ion 

After  Tr 

an:;  f 

urmat i on 

- .05702 

+ 

10. 34654 

-  .  0  7953 

+ 

10 . 20562 

- .0701  I 

+ 

13.96812 

-  .0698  3 

+ 

13.96671 

- .06244 

+ 

10.94332 

-  .  049  54 

+ 

10.93907 

- .04388 

+ 

8.99752 

-  .05750 

+ 

8 . 96044 

-1.51751 

+ 

5 . 50653 

-1.61259 

+ 

4 . 90962 

-1.61639 

+ 

5.41038 

-1  .  22750 

+ 

4.97987 

- . 27191 

+ 

1.20414 

-1  .  54666 

+ 

5  .  50147 

-  .43310 

■f 

1 .09745 

-  1  .  58461 

+ 

5 . 40713 

- .83306 

+ 

3.65462 

-1 .45632 

+ 

5 . 53310 

-1 . 22213 

+ 

3.06310 

-1.51763 

+ 

4.45319 

-1.42714 

+ 

5  .  58897 

-  .  54310 

+ 

3 .63429 

-1 . 58679 

+ 

5 .40001 

-1  .  06302 

+ 

3 . 29105 

-  .71819 

+ 

1  .2  7304 

- . 36782 

+ 

1  .131.01 

-.36226 

+ 

1.51909 

-  .  26965 

+ 

1  .  18851 

-1.06186 

+ 

2.50015 

- .43999 

+ 

3  .  85085 

- . 73434 

4- 

3.05232 

-1.25745 

+ 

3.65065 

-1 .43137 

+ 

3  .  35652 

- .62836 

+ 

1  .  31448 

-1.08042 

+ 

3.88162 

- . 34492 

-t- 

1  .  51727 

-1.61259 

+ 

4.90962 

- . 73034 

+ 

3.04703 

-1.61266 

+ 

4.90964 

-1.09035 

+ 

2 .51427 

51 


C  rntrol ler  // 1 
1 ,4 ,6 , 7 


Re  s  i  dna  1 
9,10,11,12 


Table  Xa 


'lime  Response  -  Two  Controllers 


Modal  Assignments 


Controller  #2 
2,3  5.8 


Before  Transformation  After  Transformation 


Time 

Los-X 

Los-Y 

T  ime 

Los  -  7. 

Los-Y 

•  5 

.003596 

.001203 

•  5 

•  003751 

.  >012  92 

1 .0 

.003360 

-  .000289 

1 .0 

. 00 5040 

. JO 07 04 

1.5 

.001491 

- .001148 

1.5 

.  0  02  6  2  0 

- . 000367 

2.0 

.OOOI63 

- .001043 

2.0 

-  .  000043 

-.001124 

2.5 

- .000598 

-.000703 

2.5 

-  .000760 

- .000740 

3.0 

-.000939 

-.000256 

3.0 

- , 001246 

- . >00362 

3-5 

-.000707 

.000278 

3*5 

- . 001 542 

- . >001 54 

4.0 

-  .000198 

.000628 

4.0 

- .000516 

. 000464 

4.5 

.000143 

.000625 

4-5 

.00041 5 

. 000806 

5.0 

.000216 

.OOO376 

5.0 

. 000402 

. JO 046 9 

5-5 

.000143 

.000054 

c.  c 

J  •  J 

. 000324 

. 000144 

6 . 0 

.OOOO35 

.000200 

6.0 

.000277 

- . 000066 

6.5 

- .OOOO39 

- .000304 

6.5 

-  .  00004-6 

- . JOO336 

7.0 

-.000050 

-.000263 

7.0 

- .000219 

-  .  >00384 

7.5 

- .000013 

- .OOOI36 

7.5 

-.000093 

-  .  000201 

8.0 

.000031 

.000004 

8 . 0 

- . 000012 

- . 00002 9 

8.5 

.000049 

.000100 

8.5 

.000017 

.000064 

9.0 

.000030 

.000126 

9.0 

.000065 

.000167 

9-5 

-  .000008 

.000095 

9.5 

.00006 7 

.  000150 

10.0 

-.000037 

.000040 

10.0 

-  .000021 

. 000061 

10.5 

-.000041 

-  .000010 

10.5 

-  .000051 

-  .  >00006 

11.0 

-  .000026 

- .000041 

11.0 

- .000042 

-  .  000044 

11.5 

- .000002 

- .000046 

11.5 

-.OOOO36 

- . 000064 

12.0 

.000021 

- .000032 

12 . 0 

-  .  000002 

-.000048 

12.5 

.000030 

- .000011 

12.5 

.000037 

-  .000013 

13.0 

.000021 

.000005 

13.0 

.OOOO36 

.000007 

13-5 

.000006 

.000013 

13.5 

.000020 

.000016 

14.0 

-.000007 

.000016 

14.0 

.  00000‘, 

.  O00019 

14.5 

- .000016 

.000011 

14.5 

-.000013 

.000012 

15.0 

- .000016 

.000005 

15.0 

-  .000024 

.000002 

15.5 

-  .000008 

.000000 

15.5 

-  .000017 

-  .000002 

16.0 

.000002 

- .000003 

16.0 

-  .  00000  5 

-  .  000003 

16.  5 

.000006 

-  .000005 

I6.3 

. 000004 

-.000004 

17.  > 

.000008 

- .000005 

17.0 

.000010 

- .O00003 

17.5 

.000O07 

- .000003 

17.5 

.000011 

- .000001 

18 . 9 

.000002 

-.000001 

18.0 

.000005 

- . 000001 

18. 5 

- .000003 

.000000 

18.5 

-  .000001 

-.000001 

19.  ) 

-.000004 

.000002 

19.0 

-  .000004 

.000000 

19.5 

-.000003 

.000003 

19.5 

-  .000003 

.000001 

2).  ) 

-.000002 

.000002 

20.0 

- .000004 

.000001 
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Cone  1  us  i  on ;; 

During,  this  study  it  was  determined  that  system 
response  can  be  greatly  reduced  through  the  implementation  of 
an  additional  controller.  As  became  evident  in  the  suppres¬ 
sion  portion  of  the  research,  the  designer  can  greatly  reduce 
the  computat iona 1  requirements  through  the  use  of  structural 
symmet  ry  and  sensor  locations.  Bv  assigning,  modes  to  be  con- 
trolltd  according  to  orthogonal  grouping  of  their  modal  ampi- 
tudes  associated  to  each  sensor  location,  the  system  will  be 
inherently  decoupled  as  earlier  explained. 

In  all  of  the  test  cases  run,  the  residual  modes  were 
not  adversely  affected  by  any  of  the  control  or  transforma¬ 
tion  techniques  applied  to  the  overall  system.  As  a  result, 
including  only  the  lower  freqency  modes  as  controlled  modes, 
has  pr  >ved  valid  for  the  modeled  system. 

The  capability  to  control  the  system  may  be  increased 
through  additional  sensors,  but  it  must  be  noted  that  the 
system  will  not  be  able  to  suppress  more  modes  than  sensors 
as  was  noted  in  the  transformation  section. 


Re ommendat  ions 

The  primary  thrust  of  this  invest i;mt ion  was  toward 
the  evaluation  of  a  system  which  implemented  tv^Acnt  ral  i  zed 
controllers.  The  results  presented  indicate  the  mathematical 
advantages  of  applying,  this  technique  to  the  model  chosen. 

The  imnortance  of  evaluating  the  entire  modal  analysis  became 
evident  through  the  analysis  of  the.  modal  amplitude  vectors. 

This  single  area  has  presented  itself  as  a  key  to  real  world 
application  of  decentralized  controllers.  The  importance  of 
the  location  of  the  sensors  and  actuators  that  are  used  to 
control  the  structural  motion  is  an  important  design  tool  in 
achieving  desired  system  response. 

The  next  logical  step  in  the  study  of  this  control 
problem  would  be  the  experimental  evaluation  of  the  techniques 
applied  in  this  study  to  determine  the  feasibility  of  the 
implementation  of  the  system  described.  This  would  include 
the  evaluation  and  determination  of  computing  capabilities 
required  to  achieve  the  results  which  have  been  put  forth  in 
this  paper. 
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NASTRAN  Anal  vs 


L 


1.370 


3.505 


Frequency  and  *  c  d  e  Sheets 


Nominal  Case 


ui2  =  1.467  ui  3  =  2.965 


-2.471E-C1 

“3.999E-CI 

~6 . 368E-02 

2.746E-C2 

4. 279E-02 

2.309E-01 

3 . 678E-02 

-4.758E-02 

1.4S1E-06 

-1.48SE-01 

4.C00E-01 

-2.249E-05 

-1.963E-02 

8.329E-C2 

1 . 984E-01 

-1.716E-C1 

3.390E-O2 

4.3C8E-02 

1.14  5E-01 

2 . 977E-01 

-7 . 213E-02 

*2  ' 

6, 813E-02 

?3  = 

2 . C1CE-01 

*4  = 

-6 . 817E-05 

-3.607E-02 

7.C9CE-02 

1.548E-01 

-2. 512E-01 

4 . 347E-02 

2. 253E-C2 

6 . 8C4E-02 

3  „  436E-01 

4.397E-02 

-4.721F-C2 

9.732E-02 

-8 .19GE-02 

-1.962E-02 

5.4S1E-C2 

1 . 363E-01 

-1 . 718E-01 

5.296E-02 

4 .336E-C2 

1.0C0E-01 

3 , 894E-01 

4.397E-02 

-4.722E-02 

9.784E-02 

—  — 

8.192E-02 

_  M 

3.848 

4i6  = 

s.iso 

“7  = 

5.676 

II 

CD 

3 

5.680 

-8.783E-02 

1.353E-05 

-2.661E-C2 

"-2.994E-02" 

-5.07CE-02 

1.218E-11 

4.6C7E-C2 

-1 . 731E-02 

-1 . 299E-01 

3.402E-11 

3 . 302E-05 

8.7S4E-G2 

3.095E-01 

-2 . 041E-01 

3 . 374E-02 

4 . C71E- 02 

1.786E-01 

2 . 535E-C1 

-S.e44E-02 

2. 36CE-C2 

-3 . 514E-C1 

*e  = 

-6.057E-C6 

4>7  = 

3.231E-05 

$  0  = 

3.554E-02 

2. 866E-01 

-2.041E-01 

2.733E-02 

2. 7421-02 

1.224E-01 

-3 , 535E-01 

-5.481E-02 

2 . 798E-02 

1.139E-02 

1 . 086E-G4 

-4.913E-01 

-4.875E-C1 

2 . 494E-01 

4 .  C82E-C1 

3 . 382E-02 

3."99E-C2 

1 . 868E-01 

6.802E-1C 

-5.108E-02 

9.61CE-C3 

1.14CE-02 

5.065E-1C 

4.908E-01 

-4 . 87SE-01 

l—  - 

—  _ 

i_ 

—  mm 

8.940 

“10  = 

10.303 

“11  " 

10.923 

“12  = 

13.966 

p 

_  — 

_ 

—  “* 

9.907E-02 

-3 . 390E-C3 

" 6.370E-02 

3. 206E-02 

5.720E-02 

5 . 85CE-03 

3 . 678E-02 

1.851E-C2 

1.729E-G. 

-1.605E-05 

3 . 588E-02 

6.438E-02 

1.076E-C1 

-2.286E-01 

-2.401E-01 

-4.C26E-C1 

6.213E-02 

3 . 360E-01 

-1.385E-01 

-2. 224E-C1 

-4.953E-01 

$10  = 

4 .964E-05 

$11  = 

-2.605E-01 

*—* 

K> 

II 

-1.305E-01 

-1.679E-01 

3.783E-01 

-8 . 626E-C2 

3 . 204E-01 

-2.198E-01 

4.554E-02 

3.9U4E-01 

-1 . 587E-C1 

-1.11CE-02 

-1.471E-n2 

5 . 97CE-03 

-9.27SE-C3 

-2.743E-01 

-2.286E-01 

2.984E-01 

2. 272E-02 

-3.554E-02 

-3.C49E-C1 

-2. 719E-01 

3 . 568E-C1 

-1.109E-G2 

1.472E-02 

5.971E-03 

-9.281E-03 

■*  — j 

—  — . 

L_ 

F r  e  c  u e  r.  c  v  ar.d  M  o  d 


“1  = 


*1  = 


us  - 


*5  = 


dig  = 


♦  9  = 


v  r-_  r  .. . 


Perturbed  Cacs 


1.342 

uj  2  ~ 

1.665 

“3  - 

2.891 

“4  = 

2.957 

-3.444E-01 

5.429E-01 

"-4 .421E-02" 

_5.726E-02 

5.964E-C1 

3.135E-01 

-2. 944E-02 

-9.915E-C2 

2.330E-06 

-1.G90E-G1 

2 . 788E-01 

-1.466E-04 

-3.107E-02 

1.263E-G1 

3  ,125E>-01 

-1.760E-C1 

5.379E-92 

7. 292E-02 

1.305E-01 

3 . 046E-01 

-1.111E-05 

92  ~ 

9.756E-02 

t  3  - 

6 . 519E-C2 

?u  = 

-7.205E-05 

-5.079E-02 

1 . 096E-01 

2 . 727E-G1 

-2 , 386E-C1 

6 . 51SE-02 

4.162E-02 

1 . 2'74  E-01 

2.409E-01 

6 . 380E-02 

-6.728E-02 

1 . 571E-01 

-9  .1 57E-02 

-3.1G1E-02 

9.092E-02 

2.466E-01 

-1.75SE-01 

7.656E-02 

7.425E-02 

1 . 726E-01 

3.771E-C1 

-6.379E-02 

-6 . 726E-02 

1 . 372E-01 

9.149E-02 

3.398 

0)6  = 

4.205 

0)7  = 

4.662 

4.755 

-1.369E-01 

”:.7C6E-05 

5.571£-0~ 

17.584E-C': 

-7.906E-02 

2.487E-11 

-9.647E-02 

-4 . 380E-02 

-3.441E-01 

6 . 986E-11 

-2.248E-G5 

1. 337E-01 

1 . 621E-C1 

-2.041E-G1 

-3.44QE-02 

4 . 701E-02 

9.355E-C2 

3 . 53EE-01 

5.96CE-02 

2.722E-C2 

-4.969E-02 

$6  = 

5 .160E-06 

9  7  = 

-2.905E-G5 

II 

CD 

-*> 

■3.7811-021 

1 . 620E-01 

-2.041E-C1 

-2 . S82E-02 

3 .  E71E-G2 

7 . 309E-02 

-3 . 535E-C1 

5.644E-02 

3.245E-02 

-7 . 571E-02 

1.G03E-04 

4.872E-01 

-4.69EE-01 

1.444E-01 

4.082E-10 

-3.447E-02 

4 . 6  55E-C2 

1. 037E-01 

7.861E-10 

5 . 318E-02 

1 . S65E-02 

-7.57GE-02 

6.086E-L0 

-4.872E-C1 

-4 .698E-01 

8.539 

0)10  = 

9.251 

“11  = 

10.285 

“12  = 

12.905 

r—  — 

r-  _ 

_ 

1.445E-C1 

-5 . 777F.-03 

1.594E-0r 

8.369E-02 

8.347E-02 

9.965E-03 

9 . 205E-C2 

4.833E-02 

2.702E-01 

-3.372E-C5 

2.  :3BCE-Ci 

1.587E-01 

2.125E-01 

-2.242E-01 

-1 . 51FE-Q1 

-4.059E-C1 

1 . 228E-01 

1  3.883E-01 

-8 . 756E-02 

-2.343E-01 

-3 . 26GE-01 

$10  = 

4 . 517E-0  5 

9l  l  = 

-3.117E-01 

$12  ' 

-1.611E-01 

-1.414E-C1 

1  3.846E-01 

-1 . 619E-01 

2.996E-01 

-3.096E-01 

'  3.6811-02 

3 . 311E-01 

-1.419E-C1 

-1 . 504E-02 

-1.184F-02 

9.133E-04 

-8. 2CCE-C3 

-3.389E-01 

-  2 . 241K-01 

2 . 057E-01 

2.6S7E-C2 

-3.22eE-02 

-2.iu'n'.-ci 

-3.058E-01 

3.304E-01 

-1.5C3E-C2 

1.135  :-02 

l-  — 

9.153E-04 

-8 . 2G3E-G3 

PROGRAM  THESIS 


A  h  /  7  • 


CPT=  . 


FT  N  5 ♦ 1 +E  26 


PROGRAM  THESIS 

•IE  AL  f  ODE  (t , it ),  XL  (•»?■>  ,IN1T  <■  ,i;) , XD(4k  ) 

REAL  FSST  ( 1 2 , 5  2)  ,-UC  ( 12,  1  2)  ,  A  03  2  ( 3  2,  12)  ,  t  <ZZ  (12,3.2)  ,  GAlNi(  -2,  . 
REAL  DT  »  E  AT  1  (4  , 4  j ) ,  SAT  (  4  ,  A .  )  ,  E  U  2  (  *  ,  ,  4  '  ) 

REAL  ACT<12, 12), TT (12,12)  ,7£7>,,H(12) 

REAL  KT2(12,i2),03T(l2,I2),CSTR(it,i2),TRT(i2,:t) 

REAL  f  T  (12, 1  2)  ,  TC  I  (1  2 , 12  )  ,C  T  7  (12, 12),  £^1(12,12) 

REAL  A:C2,l2)  ,1(12)  ,FH15  C2,12»  ,::(i2,!2) 

REAL  CAT(iZ,*2),;3G(l2,J.2),S^l(12,L2)  ,  ERG  ( 1 2 , 1  2 )  ,  OR  ( 1 2 , 12  ) 
REAL  KA  JR  <-.,<•  )  ,GAIN(12,l2)  ,S:3C2,lt)  ,A*C(i2,it) 

RE  AL  KTS(3t ,32) ,33G(l2, 12 ) ,Co6(lt,  t2)  ,  CT  (  12,  it ) 

REAL  V(it,*2)  ,  S±t*G(iZ) ,  I  R  C2,i2> 

REAL  C7C:  (1 2,1 2) , P03 (12, 1 1) , ACG( 12,  12) ,<T1 (1 1, : 2) 

REA  L  XTR(l2,l  2)  ,  310k(12 , 12)  ,11.  ,  F£\!(  12, 12)  ,  AA  ,?E,:A  ( It,  12) 
REAL  rsi  (12,12)  ,7  1(12,12)  ,7  7  (12,12)  ,R  (It,  12)  ,ki  (12,12) 

REAL  c:  (12,12)  ,phi  (3  2, 12) 

INTEGER  K,uC  ,NS,N  7,IC  (12)  ,i  C2,NR2,  (52  ,MM,L,P  ,  DtC,  E 

INTEGER  I,)  Erv,J,'TZ,.'i2,0,I  *>PE  ,  13(12), N ACT,*)  (12) 

COMPLEX  H 1  ( 1 2 )  ,7(  ,.) 

REAL  K0£(i2,12)  ,  KCC  ( 12, 1 1 )  ,  WL  hm  1 , 4 ,  > 

COMMC!  /t  t  !*■  i/?,DIi',NOlMi,  T  El. 

COMMCN/hA  i M  /NCm,  NOAl,WjFK 
CD  MM  J).’/)' A  IS  2/3  TOc 
CCMMON/f  Ai  KJ/XT? 

COMMCT/INOOT/T  APf 
COMMOf./N'J.H/IC,  JS,IR,NC', \S  ,NR  . 

CCMMON/SAVi/7  (1J  .)  ,TS(i.  k  ) 

NDIM=12 
N  0 1 M  1=1  3 
NO  A=  r 
NO  Al  =  33 

n= : 

TA  CE  =  9 

PRINT*  ,  'EUTcn.  t  C,NS,  Nf.,N*C7  ,  US:  H,  ?i.T  A>  * 

RE  AD  ■  ,NC ,  NS ,  17s  ,  NACT, NSEN ,  7:K- 
PRINT*,*  * 

PRINT*,*  EnIlR  THE  *,NA:r,*  lLEHENTS  FOR  EACH  -MIA* 

_N  =  ND  +I-?  +HR 
DO  S  7=1, N 

PRINV  ,  »£.-,TEr.  c H 1 A  *,I,»>* 

\cAO(G,’)(FHl  (I, J) , J=l,*,C7) 

9  CD  NT  J'.  IE 

=  RI  JT*  ,  '  * 

PR1N'*,*  _m1.:  the  *  ,  nil-E  N ,  *  tL-HENfS  FCR  E-CH  -HIS* 

DD  6  1=1,  N 

PRINT*  ,  *  E  Mi  E*  PHI  3  *,I,*>* 

REA  IC  ,  •  )  (FH1E  (i,  J)  ,  J=mMLi.) 
fc  CONTINUE 

QD  h  1=1, N 

PRINT*,*  ENT  wk  OMEGA  *,[,»  >* 

READ  (c  ,  *  )  W(  I ) 

0(  I)  =  -2*  ?E  T  A  *  V.'(I) 

4  GCNT  :?M!E 

c 

I, 

C  INITIAL  Duiif  tTIJ.NS  AkE  ncAL  ih  3 f  RDW 

C  ROW  1  IS  X  ,  •  UW  2  IS  X  j  CT  ,  \D  H  J  lb  E,  -OW  IS  E  DOT 
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PRCGRAH  7  He  SIS 


7h/7* 


CPT  = 


PTN 


5.1  +  S 


C 

c 

00  14  1=1,2 

19  REAO(f,M  (VC0£(I,J),J  =  1,N) 

00  21  r  =  l,<; 

21  R£A0(£,*)  (IMT  (I,  J),  J  =  1,N) 

299  CONTINUE 
0£C=  - 

PRINT*,*  IF  7IS  IS  A  OEClUF-  LLD  RJN  ENTER  1  ELSE  ENTER 
RE  AO  , OEC 

PRINTS*  0'_ COUPLE  =  *,0-0 
IF(O.EO.c)  7  WET 

PRINT*,*  ENTER  NC,NS,nF  >  *. 

RFAO-,NC,NS,NP 
EM  Cl  c 

PRIN7*,*  ENTER  THE  *,NC,  *  Cul;T  ROLLED  NODES  >  * 

KE  AO ’-,(10(1)  ,1=1,  NO 
PRINT*,*  *  ,  CG(I)  ,1  =  1, NO) 

PRINT*, 'ENTER  THE  *,NS,'  SUPRlSSEG  HODES  >  * 

READ-, ( IS (i) ,1=1, NS) 

PRINT*,'  *, (13(1) ,1  =  1, Ni) 

IF  (iiF.NE.  )  THEN 

PRINT*,*  ENT  c.K  THE  ' , NK,  '  rSSILJA.  HOCES  >  * 

RE  AD  S(Kd)  ,1-1,  NR) 

PRINT*,*  *  ,  ( I  <  ( I )  ,  I  =  1 ',  T  R ) 

en  r.i c 

NC  2=  2*  NC 

NR  2=  2* NR 

NS  2=  NS*  2 

N2 =2  'N 

CALL  FOR'iXO  (XO  »1MT) 

IF  (DEC.  EC.;)  CALL  cC-HXi  (X  0,1  ill ) 

M=  2’  :IC2  +NS1 +NR2 

IF  (DEC.  EC.l  )  +  c+t  Kl 

print*,*  Initial  conch  k  ns* 

CALL  PRNT(XC,M,i) 

PRINT*,*  TO  FRi.NT  ALL  uF  THE  HMR101ES  E'iTEF  1,  ELSE 
REAP  ,n 

IF  (0.“0.1)  THEN 

p r i n 7 * « .*  t h {  a  control  'iatox  is* 

CALL  F  P  F  n  A  ( l  0,  r,W,iO  ,NC2,  1C) 

CALL  Fn M  (AL,NC2,.,'C2) 

PRINT*,  «  -HE  P  CCMROLLEL  IS  * 

CALL  Ff)F  (•  8  (  Ec  ,  F  HI ,  i,C  ,  NC  j: ,  HA  C  7  , 1 0 ) 

CALL  F-  NT  (EC  ,  *  C  2,  4'  ACT  ) 

PRINTS'  THE  C  CONTROL. EL  R,7*IX  IS* 

CALL  PORnC (CC,PHIS,NC,\C2,NSlN,IC) 

CALL  ASM  (1C,  *.S_N,NC2) 

PRINT*,*  THr.  I  SUPRESSED  t  A  7  r  1  <  IS* 

CALL  FOR^A  (A  C,  0,W,u5,NSE,  IS) 

CALL  F-N!7  (AC,RS2,NS2) 

PRINT*,*  7HL  E  STRESSED  f i  r.  1 X  IS* 

CA  LL  PORTS (  EC  ,  PHI,  r.3,  NCI ,  f:»  C7  ,IS) 

CALL  f  1.1  (jC,.O0,  .017) 
p  ••'•,*  M  .  1  •.  :s* 

CALL  Frrvi  C((  C,  F  H 1 3 ,  N  S ,  t.  S  2  ,  !•  6  .  h ,  I  5 ) 

CALL  Pf  M(CC,l.i_N,NS2) 


P*1N7*,  '  7  He  c  HAT.iX  15' 

I  CALL  FORN3(GC,rH:,UR,NRc,NAO-j  ,1R> 

1  CALL  FRNT (3C, N? 2, NACT) 

PRINT*,'  iHi  C  RESIDUA  _  -Air  IX  13' 

:  CALL  FORr  C(CC,  FHI3 ,Nk, NR£  ,f  SEN,IR> 

*  CALL  PKNT  (CC,  NSEN,NR2) 

‘  0=  t- 

ENDIr 

i  CALL  FOM  3(cC,  PHI,  NC,N,'2,h;Cl  ,J5) 

!  CALL  TPR  (B7 ,3C  ,:JC2,MACT,  l  ,2) 

•  call  r*u»  ("c, z:  ,NC2,r<Acr  ,nCc ,sar ) 

1  CALL  FOR)  C(CC,  FHIS,NCf?.C2,t.LlN,i:) 

I  CALL  iPR(Cl,Cv,NEeN,NCc,A,£) 

;  CALL  f'^UL  (CT,  CC  ,N02,  NSLN,  f.C£  ,  CTCO 

i  12*.  CONTIMtE 

\  72 

PRINT',*  ENTER  THE  DIAGONAL  TER  4  "DR  THE  HEIGHTENS  MATRIX  C.  >' 

RE  AD  S A  A 

printsaa 

PRINTS'  ENTER  THE  OBSERVER  W  Ei  GHH  NS  DIAGONAL  TERN  >' 

RE  AO  “»PF 

PRINTS  EE  ‘ 

DO  1C  l  1=1,  NG2 
DO  i (  J=1,NC2 
IF(I.FSJ)  THEN 


IS 


DA  (I, J) =:A 
QO  9(1, J) s 3? 

ELSE 

DA  (I,  J)  = 

QD?(  T.,J)  = 

EM  Cl  ~ 

CONTINUE  < 

CONT  I  MJL 
IE  c= 

TO  L=  •  i 

CALL  FOR*' A  ( A  C  ,  D,  W ,  MC ,  NC2  ,  It) 

CALL  F  (NT  2,  ,,C,SAT  ,  Om,  LaT  ,  ^  ES,  I  0.,  ICR) 
j. F  rs SO.  )  T  HeU 

FRIN“S'  THE  \ I  GATT  I  SOLUTxLt.  G r  AS  PCS' 
PRINTS'  IlF.s  SIER 
CALL  r»M  (CAT  ,f  C:,NC2) 

ENCIP 

IEP= 

TDL=.  1 

CALL  7'*T(«CT,AC,NC2,UC2,1,2) 

CALL  F  r- 1 C  ( NC  2 ,  C7  ,  CT  CC ,  0  0  ° ,  Fo  &,  A  3  3 »  T  D  L  ,  I E  R) 
CALL  MMUL  (CC,FC3,  ^SEN,NC2  ,NS2,Kf  i) 

IF  (7  -J,  -Q.  1)  7 m E N 

CA  LL  HMUL  (■■  1  i,‘,  ,RT  ,F,  F,N;Et.,ST  OR) 

CA  LL  F  *'L'L  ( >  T  G'  ,  RT  i,  3  ,  ,<S^  M  ,*  C  u  ,  <F  2  > 

CA  LL  1,'UlU(  ,<T  2,  S  \G_,  TT  i) 

CALL  7 rF  (KCGjKT  1,NL£N,N:2-,:  ,  2) 

PR  IN?*’ ,  '  THE  <  GAIN  Ntit.X  ' 

64  •  • 

■•S*  » 


I 


r 


T-iiSIS 


T4/7V  LPT  =  . 


F7S  5.1  *ZZi 


CALL  FK  NT  (KC3,  r.C2,  NS9N) 

CALL  F“'Ul  (Kt3,C3,Su2,N<?LNff.CcfIC3) 
CALI  FOK^CtCC,  f4i  3,Ni>,K5  2  5) 

CALL  K'lUL  (<C9,CC,  IL2  ,  NS  -  ri  ,f,b  2  ,  *3  5  > 
CALL  FORmC  (  CC  ,  FHI3 ,  NF?  ,N  \  2  ,1  S  -  K ,  I  ^ ) 
CALL  MA-JL(^OF,CC,NC2,NS£L,Nk2,<:^» 
QO  S’7  1  =  i »KC  2 
00  5?  J  =  1 ,  N  C  2 

AKCC,J)  =  AC  <  2 , J )  -  KCC(I,J) 

HM=NC2+NC2  +  f.S2+‘J^2 

xF  (09C.  f0.i)NM=2'  NC2  +  2-  NS2+M\2 

00  9  1  I  =  1 ,  h  A 

00  9 i  J=l,Mh 

HA  JM  (I,  J)  =  . 

OO  92  I  =  1 »  H  C  2 
DO  92  J=i,NC2 
MAJM  (1,  J)  =;  CGC>,  J) 

CA  LL  FOF  f  -3  (  FC  ,  PHI ,  NC  ,  NC2  ,  N/  L  7  ,  J  0) 
CALL  TFF  {  9T  ,3C ,NC2,NACT  ,1,2) 

CALL  KMUL  <P7  ,  C  A  T  ,‘IACT  ,A.X  2  ,AX  2  ,  GU  M ) 
DO  T=i  ,NACT 
DO  77  J=1,NC2 
GAIN  (I,  J)  =-G^I  A;  (I ,  J) 

CALL  >MUl  (FC,  C-AI N  ,  NC  2,  M  A-CT , i  C2,  32  2) 
CALL  FOF>  9  C  CC ,  FHI ,  A<3  ,  NS 2  N.,L I  ,13) 
CALL  *MUl  (vC  ,  G^I  J ,  N  O  2  ,  M  C  7 ,  N  i.  2 ,  ^  3  3 ) 
CA  LL  FOP  f.3  (Cl  ,  f  Hi  ,  u  <,  N-.2  ,  l.-'  CV,* \) 

05  LL  MwU-  (r,C  ,G;I..,Na2,SAC7  ,KC2, 3\3> 
L=2‘NC2 
no  93  1=1, NC2 
00  93  J  =  I ,  C  2 
iA JH(I, (J+„C2) )=  3CG(I,J) 

30  3  i  I  =  i , A  C2 
00  9-  J=  1,  AC2 

tlAJK(Tf\C2)f  (  J+NC2)  )  =  AKc(l,J) 

00  9:‘  I=l,i.C2 
00  9;-  J=i,N5  2 

M A  JH  ( ( I  +  f:C2 )  ,  (  J+L  )  )  =KC3(I,J) 

00  93  I=i,N$2 
00  9^  J=1 , NC 2 

M A  J’l  (  (L  +  j  )  ,  J  )  =  SSG  (  i  ,  J) 

!A  JH  (  (L+I) ,  (Jt-?C2))  =  9SG(i,J) 

CALL  FOF  A  A(  AC,J,W,NS,NS2,15) 

DO  97  1=1, M2 
DO  3”  J  =  1 ,  N  £  2 
M6 JM( ( I +L ) , (J+L))=AC(1,J) 

.1=  L+NS2 

CA  LL  F OK1'  A  (/  L ,  C ,  W ,  Nrt  ,  Nk2  ,  ) 

DO  3  1  =  1, Nf  2 

00  3  J  =  1,NF2 

M6  J'«(  (1+-  )  ,  (  J  +  M  )  =  AC  (I  ,J> 

00  31  l  —  1 , A'* A  £ 

DO  3  i.  J=l,‘JC2 

:A  J.H<  (I  +  v  j  ,  j  )  rp,_&<  jj 

■ f  J"  ((:♦'),  u+-  32) )  i,  j) 

'  A  JM  (  ( J  +  \C2)  ,  (.!+*■))  =  <Z^U>.) 

IF  (0  ;C»  £.0.1)  T  h: £M 


1 ,  N  C  2 

o  =;  cgc  ,  j) 

A-3(  EX,  PHI ,  NC,  NC2,  W  C7  ,JC) 
(9T ,3C,NC2,NACT,i,:) 

L  <P7,CAT,NACT,A!C2  ,AX2,GU 
1  ,  N  A  C  T 


ac  a  ,j> 


=  EV.&  (  *  ,  J  ) 

= "  ■  r  c  ^ .  j ) 

—  \  (  j,  * ) 
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th£s:. 


7h/ 7, 


CPI  =  \ 


FT S  ».l  +  :26 


DO  4  >  1  =  1 , NR 2 

OD  4  J=i,4C2 

MA  JM  (  (I+M),  (J+EC2)  )  =-.: 

M4JM  (  (I4M)  ,  J)  = 

U  MA  JM  (  (J+MC2)  )  =  w... 

CALL  FOrrA(*.C,  C,W,NS,NS2,  II') 

CA  LL  FOM  B  ( EG ,  FHI ,  1.3  »  N52  ,  Ut  C*i  ,1  5) 

C4  LL  7rM  3T  ,5-  ,.4E2,NACr,  i  ,c> 

IF  (77«:Q.J  CALL  MMJL  (£3  ,  £1  ,  i\S<L,  Ml  3  T  ,  NS  2 ,  3S6T ) 
IE  R= 

TOL=.:  1 

CALL  I’RiC  ( KE2»  AC »  TSTT  » CA  ,MC,A£-:S2,rD-,IE?> 

CALL  TFK(ACT,AC,'JS2,NS2,1,2) 

CALL  FCiRKCtCC,  FHIS,NS,NE2  ,;<S"NtI3> 

CALL  TFKCi  ,CC,Ni£N|NS2,i,2) 

CALL  M*UL <CT,CC,N32,NSEn,Mc,CTCC> 

CALL  K.'i:  (uS2,/CT,CTCC,2C£,FCE,AC3,  T3L,I:R) 
CALL  Tck(AKC2, ACG, NS2,KS 2*i, 2) 

M® 2*NC2 
DO  41  1=1, NS 2 
DD  41  J=1,NS2 

MA  JM  (  <M+i  )  ,  <M+  J) )  =  A3G2(I,J) 

41  MA  JM(  («  +  \S2+i)  ,  (V+NS2  +  J)  )  =  ;KC2([,J) 

CALL  ►•UUL  (CC  ,PC0,  JS:.N,NS2,I  G2,Kri) 

CALL  TFRtKOE  ,<Ti,  JSEN  ,!o  2  ,  2) 

C  KOB  Ij  ‘IOW  THt  •<"  GAIN  MATrwX  FGF  S  f  3  T  E  M  2 
M=  H+.4S2 

CALL  FQR4GCCC,  FHIS, HC, NCt ,4$ s K, I 3> 

CALL  M"UL (KC8,CC,NS2,NSiK,NCL,<:C) 

DO  42  1  =  1, NS2 
DO  42  J  =  l,l,C2 

i-  2  MA  JM((M+1),  J)  =  *<CC  ( I ,  J) 

CALL  MU.  (  0  T  ,  n  I C ,  ^GT  ,  N5  2 ,4S  2  ,  G  A I  S  2) 

DO  7.'  I  =  1 , »;A. C 1 
DO  7  3  J  =  1 , NS  2 

7c  GAIN2(I  ,  J)  =  -Gf  INC  (I,  J) 

IF  <7  7. f.Q.  1)1  Hi 4 

CALL  M'-U_(T*i  ,G*If-2,E,NACT  , t  S  2  » S  T  0  F. ) 

CALL  MMUt.  (^l,*>Tj  ;,£,:.,N52,Ui.) 

CALL  MU.  (Tl  ,  i  i.  N ,  N  A  C  T  ‘it.Z , A1  >J  2 ) 

ENCI- 

CALL  r0M  5  ILC,  FHI  ,NC,MC;  ,  ML'  ,12) 

CALL  7  MU.  (8C  ,GfI,;2,.)C2,  'U  Cl  ,i.EZ,  303) 

M=  710  2*2 
DO  43  T  =1 , 4C  2 
00  4  3  J=1,NS2 

ma  jm  a,  ( ^  +  j) )  =  ncc ( i , j; 

43  MA  „M  (J,  <*:♦»,«  2*  J)  )  =  £CC(1,J) 

CALL  F03  *  3  (EC,  -HI ,  i<3  ,  N?  2  ,  N;  Cl  ,13) 

CALL  MMUL  <6C,GAIN2,NS2,SAcT,r>S2,3SG) 
j  0  ►  *  I  =  1 ,  N  f  c. 

DO  4>  J  =  1  , f»  1  2 

44  MA  JM(  (M+]  > ,  ( r-ST  J>  )  =  3SGC,J) 

■<=  2*  If  2  +  2-  :,i'Z 

CC  -  >  T  =  1  ,  \ F 2 
.  0  '  J  =  *  i  <'C  c 

MA  JM<  (M  +  I)  ,  J)  =  ?-,G(I,J) 
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i  PROGRAM  H; SIS 


74/74 


CP7  =  i 


PIN  5.1-C2* 


V  ^  k-  t  r  ■  I\i  '  w  i  W  ;  f 

CALL  MMUl  (KC'SC 
00  4->  1  =  1, KF  2 
DO  4 i  J  =  1 ,  N  S  2 

Li  A  III  1  <  LI  .  «  \  •  I  i  ^ 


DO  4 i  J  =  1 ,  N  S  2 

MA  JM(  (M+i),  (  J*2*N02>  )  =  3r.CC  ,  J) 

MA  JM  (  (MS  ) ,  ( J+2*NC2*NS2)  )  =  2. KG  <  I ,  J  > 
HA  JM  ( <2*  NC2+NS2+ J)  ,  (M+I)  )  =  KCa(J,i) 
CALL  F'U-vA<  AC,  C,w,N?,Nk2,  IF) 

03  47  I =1 , NP2 

DO  4/  J  =  i ,  N 4  2 


DO  C  J  =  i , N 4  2 

MA  JM(  (M+I) ,  (M+  J) )  =  AC  ( I ,  J) 

ENDIF 

IF  (DEC.  £0.1)  MM  =  2*.NC2*2*NS2+NR2 


FORMS  £  10  THE  AT 


IF  (77.LD.1)  THIN 

PRINT"1  , *  1  HI  SUPRESSEO  ANALYSIS  [3  CALCULATE?* 

EN  01  r 

PRINTS'  F  OR  THE  TIME  RlEPLNS*  A  Y3  LIGfNVALU-  ANALYSIS  .  M  i  * 
PRINT*,'  FOR  ONLY  TH.  tIGe.f-V4.JL  ANALYSIS  2  >* 

MEADS  SKIP 
IF  (3KIP.fc0.2)  THEN 
GOTO  21- 
EN  Cl c 
OT  =  f.ll 
TOL=.rs 
03  I=:,K.M 

03  9  ■■  J=l,.- F 
EA  T1  <3,  J) =  ., 

£AT1(I,I)=1.  ? 

VO  rK  ( J,  J)  =  ‘*A  JM  (I ,  J)  ‘OT 
£AT(I,J)=WOf <(I, J)  ' 

M=  1 

LL  =1 

CONTIH't 
00  iis  1-1,  MM 
00  11'  J=1,»M 

EAT2  (1,  J>  =  EATl  (I ,  J)  +  WOi  <  (1,  J)/L. 

M=M+1  u 

LL'-LL'M 

DO  111  1=1, M 
DO  111  J= 1, f  M 

IF  (A.  3S(eiT2C,  J)-£mT1  (I ,  J  )  ) ,  G1 ,  f  0. )  THEN 
DO  112  L  = 1 , F  M 
DC  lie  K  =  i , I M 
r"  A  T1  (L,  <)=a  T2  (L,  <) 

1AT2  (l  ,  K)  =  Hr  K(L,<) 

D  “  l  L  V”  1  ’L  F-  ( :_ »  i  ?  ,  -  I  ,  •  1.  *  ,  ‘  S  k  ,  ’  ,  •i.j'fCrK,-.  ,  i  L  ) 

G J  TO  1 1 c 

:ndic 

- - *  '  -  •  -  67  . •-  . T 


» 


PROGRAM  THESIS 


7A/ 7A  C  PT  =  w 


FTN 


111 

C 

C 

C 

C 

c 

c 


21L 

C 

C 

C 

c 

c 


9- 


to 


t? 
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1 


CONTINUE 


THE  SOLUTION  TO  f.  TO  THE  AT  IN  SAT2 

THIS  FLCCK  DETERMINES  HE  LlS  A^O  3RIN'TS  THIS  VAl J 


IF  (DEC*  EC*  ■.  )  THEN 
CALL  FOPNXC  (Xi.  ,INIT) 

ELSE 

CALL  FOFrXI(XC,IMT) 

EN  OIF 

CALL  TIKE (EAT  2 ,MM,DT,X1,XG,NLDE, EA  I- ,  «f  0PK ,  D  _C  ) 

CONTINUE 


EIGEN'  VALUE  ANALYSIS  SECTION 
CALL  EIG'  F  (  Kw  JH,  rir1,4  -  » T »  Z  ,r  c  l», SOI  *1.  WORK,  Ic.R) 


POINTS1  C  V 1  f  ALL  SYSUM  tlGtN  VA.JES* 

PRINT»  ,  *  It\  =  MER 
DO  ft  •>  .1  =  1,  HP 
PRINT*,*  *,7(I> 

PRINT1,*  * 

CALL  EIGpF(ABG,UC2,NDIM,.  ,Wi  ,  IE  V  ,  vi  01 1 ,  STOR  ,  I  ER) 
PRINT*,*  EIGENVALUES  0r  AC  +  SOS' 

PRINT* ,  *  I  HR  =  * ,1 ER 
GO  9  1=1, NC2 

PRINT*,*  *,W1(1) 

PR  IN'* ,  *  EIGENVALUES  Op  AC  -  KCI* 

CALL  5IGRF{m<C,NC2,.4DI?’,.  ,«1  ,  Tc.N  ,  M  01 1  ,  STO  R  ,  I  ER) 
PRINT*,*  IS*  =  * ,  I  EK 
0  0  It  1  =  1,  NC  2 
PRINT*,*  * ,  W  i  ( I J 

IF  (0  EC.  £3.1)  THEN 

PRINT*,*  EIGENVALUES  OF  A+  EG  SYSTEM  2* 

CALL  EIGRF  ( A  TGI,  :!>2,W0Ifc. ,  „  ,  W  i  ,  T  E  N ,  T  DI  r, ,  S  T  OF  ,  IE  ) 
PRINT*,*  IcR=  *,I  Eiv 
00  6/  1=1 , NS2 
PRINT*,*  *,W1(I) 

I E 

PRINT*,*  EIGENVALUES  OF  „  -  KC  3  Y  S  T  LN  ?* 

CALL  F:GRF(AK:2,NS2,nOH,  ,Wi,T£.N,  N  0 1  N «  A,  IS") 

°RINT*-,*  IE*-  * ,  1 E  K 
DO  O'  1  =  1,  VS  2 
PRINT*,*  *,Wi(I) 

EN  01 " 

IF  <77.£Q. 1) GOT  C  2: 

CAL  l  F0r*’C(CC,(-HU,NS,NS2,f.SEN,r?) 

0  A  LL  i  F R  ( G  i  T  ,C  L  ,  Go  c. N ,  NS  c  ,  1 ,  c  ) 

00  1  f  1=1,  AS 
00  1  J=1,N"~U 

7(1,  n  =  C-E 7  (I  ,  J) 

■J 1  L  L  l  V  j  f  (  J  ,  ‘  i  x  ;  ,  i 1  <  j  ,  m  3  i  <  ,  c  f  i ,  N  j  I  •vi ,  ”  l ,  i  1 ,  >  G  ,  S  %. , 
°R1NT*  ,  *  liVTF  IEl=  *,IET 


PROGRA.-  TMESIS 


7L/7i 


GPT  = 


F7h  ;• 1+i 2  c 


IF  (F.LT.i)  THEN 
00  i  1  1=1,FSLN 
Id  T$  (I,1)=V(I,N5EN) 

P=1 

ELSE 

00  HE  I=i,NSEN 

00  109  Jrl,p 

199  TR  (I ,  J)  =  V  (1 ,  ( J  +N> ) ) 

EN  DI  F 

POINTS'  1  r  £  NSFOkMATION  J-.t  T  Fi  <  • 

CALL  FRN1  <  V  F  ,N  SEN  ,P) 

CALL  uHUL  (CSl  ,';R,NS2,NSLN,F,CSTR> 

PRINTS  *  CST  *  T  l> 

CALL  FF;NT  CETR,NS2,°) 

PRINT*,*  .'Ho  SINGULAR  VALUES  * 

CALL  Fc  NT  (SING  ,NS,  1> 

CA  LL  T  Frc  (1  :<1  ,  TP,  NScN  ,  P ,  i  ,  2) 

CALL  MUL  ( H  V ,  TR,  ? ,  fiSEN ,  P  ,f  i  ) 

CALL  C-T1KV  (  F  , p  ,R7  ,  RTl ,  J,  T  AP o ) 

CALL  F^fJ  C(CC,  FHT  $  ,NC  ,  NO  2  ,M  EN,  10) 

CALL  IFF (C7 , CC , NS IN, NC2, : 

CALL  HMUL  (TFcl  >  CC  ,  3 , N SEN,  S CC,  1  CO 
CALL  ,^!Ul  (CT  ,T  F,;,C2,  N  SEN,  P ,  C  "•  r ) 

CALL  M‘?UL  (C*.T  ,Hl,N02,F,3',Si  UK) 

CALL  FMUL  (ST  Ol  ,  TCC  ,  NC2  ,F,f.CL,CIOO) 

77=1 

IF  (0  EC.  t/3.1  )  THEN 

CALL  FOf'h  ,  UC  ,  NC2  ,  N..CT  ,iO> 

00  7  1=1, ‘,C 

00  7  J=1,NACT 
7<  V(I,  J)  =  8C  (  ( I  +  NC )  ,J) 

CALL  LSVjF(V,  !,t.C,N?iN,TL.,,NOlr,  -  l  ,  SING,  ST  C-  ,  *L\) 
PRINT*,*  LSVD"  Fui  CCxTfLL  S=lL.  OVER  I*-.  =  »,ifR 
2  =  NAFF  -  SC 
IF(t.li.l)  THbi; 

OO  71  I  =  1 , s A C T 

71  11(1,1)  =  V  ( j. ,  FACT  ) 

ELSE 

DO  i  2  1  =  1 , NIC  . 

00  7  2  J  =  l,-1 

72  T1 (I , J)  =  V (1 ,  (J  +  NC)  ) 
i  N  HI  F 

PRINT*,*  S  YSTEN  2  Tk&  NS  FO.  L  0  N  MATPI<* 

PRINT*  ,  *  • 

CALL  FRNT  (  T  2  ,  N  ACT ,  E) 

PRINT-,*  * 

PRIN“S *  .HE  SINGULAR  VALUES  Cr  =1* 

CALL  FRNT (SI N&  ,NC , 1) 

CALL  Ter  (J!  ,T  1  „j£  OT,  E,1 ,  2) 

CALL  K«U.<n,i  l,E,NACTfi,r) 

CA  LL  C-PI  (  V  (  i  ,  t  ,  R,  <1,  J  ,T  A  P  £) 

CA  LL  Fr;.‘  p(  tr,  cht  ,,N0,NS2  ,  I  ;  CT  ,13) 

Cf.  LL  *  ’Ut  (EC  ,1  ,  t,  =  :  T  ) 

u  1  LL  1  L  u  (  b  •  t  i  .» ,  ,  „  i  ) 

CALL  R*’UL  (  21 P  T  ,  IT  ,  N3  2 ,  E,  N  M  I  ,EST) 

CALL  TFR  (STCK, pC,NS2,NA0  7  ,., 2) 
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'  CALL  HMUL (PST, STOR , NS2 ,  N  A  CT  ,  NS2 , 35 3  T  ) 

IN  01 F 

I  GOTO  r 

Zt.  CONTIFUE 

'  PRINT*  ,  *  THE  CHECK  FOR  L I N  t.  A  R  G  0  1 30NA T I  0 f.S  ON  CONTROL  * 

L  =  NS  +1 

DO  3  i  1=1, NC 
.  H=IC(1> 

,  CALL  FGK.FC  ( CC ,  PHIS  , NS  ,  N5  2  ,t<S  E  N ,  I  S ) 

00  3  1  K=  i,tS 

00  3  1  t=l»rSEK 

j  3*.  1  V  (  E,  K)  =  CC  ( E , *) 

DO  3  2  J= 1 , NSEC 

3,.  2  V  (  J,L)  =  F-'IbC  ,  J) 

CA  LL  LS  VC F  (  V ,  NC'iM,  NS c.N,  L ,  T  LN  »  NL-L  1 ,  - 1 ,  SING,  Cl  C  R ,  I  c.  R ) 

'  P^INf*,*  THE  C0N10CL  NODES  USED  IN  THE  CHtCK* 

PRIN'S  (I  b(  K)  ,  H=.  ,NS)  ,k' 

PRINTS'  SINGULAR  VALUES.' 

3CC  CALL  FRNT(STnG,l,1) 

IF  ( DEC «  E  0«  3  )  THEN 
L  =  flC  +1 
DO  3  3  I  =  i,r.S 
*  M=  IS  ( I) 

1  CALL  *FCFH3(cG,FHI,NC,NCZ,NACT  ,10) 

]  DO  3  h  K=1,NC 

!  DO  3  A  E=l, r ACT 

l  3>  A  V  (  K,  S)  =  3C  (  ( K  +NC )  ,  £) 

DO  3  r  J=l,  NACT 

3-  l  V  (  L ,  I)  =  F  HI  (SJ> 

CA  LL  ISVCF<V,HCISL,NACT  ,  f  t). ,  Hut  H ,  - 1 ,  SI  MG ,  S  T  CP  , . \  ) 
i  PR  IN"* ,  *  'HE  CONTROL  NODES  USED  i  N  i  HE  CHECK' 

■  PC  INr*  ,  (IC  (K)  ,  K=i ,  NO  ,M 

PRINT*,*  THE  SINGULA'  VALJ.S* 
i  3.3  CALL  DENT  (SING  ,L,  i) 

|  E  n  n  j [r 

.  PR  IN  i‘  ,  *  TO  CHANGE  TH-E  WEIGHTING  MATRjS  ENTE'  1  * 

1  PRINTS*  TO  KEA.;nA  f.C-E  “OLcS  FOR  n  =  *,N,  •  Of.’  TO  HAKE  A 

I  PRINTS'  DlCCU-LcD  RUN  S !  T  _  k  l  * 

PRINTS'  TO  TE  r  MI  NATE  THl  )  U„  tSTER  3  >  ' 

!  READ  ,0 

PRINT* ,0 

,  IF  (T.EC.l)  THEN 

'  GOTO  12 

EL  Sc I F  (C  ,  i C.c)  THEN 
GOTO  <09 
E  N  DI  r 
END 


MAP — (LC=A  ) 


DRES3--SCC<- 

. SIZE 

-N  A*‘ : 

—  hog-.e 

>31  9 

REAL 

ACT 

c  ;  t 

-S3 

REA  . 

L  +  + 

A  W'C 

31.  <• 

,s*  4  r 

%  »*;  _ 

*  /  V 

a  hc  : 

v.  +  . 

2  ’>  i' !.  '* 

REA  l. 

Lt4 

98 

3  mu  2 

-ifct  3 

REA. 

PC 

2  t>3  f: 

"  ‘ : 
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* 

' 

> 


'4/7, 


CPT  = 


FTn  5» 1+E2E 


JSROUTINE  1 1 n£ 


SUPkOUTI  f«L  U.-  E(-  AT2,MR,OT,Xl,XD,*10DE,t4f  ,  WOFK^'tC} 
COMMjMNAINA/NCA,  ADAl 
COMMOF/StVr/T  (!■■..)  ,TS(i:  •  ) 

C0MM3f/NUM/:C(l2),:S(i2)  ,Ir  (i2),M:,^3,fn 

REAL  XO  <NL£.  )  ,£  *-T2  (NO  A  ,NL«  )  ,  E. ;  T  H  LA,  NO  A  )  ,  Of ,  MCD:  ( 2,  12  > 

REAL  WuF.K  (NCA,  F  OA)  ,  A  ,AA,  7  ,Xi  (N3A) 

INTEGER  V'M,  CEO  jZ"* 

N=  1 
«=f- 
A=  1 

2‘.  3  CO  NT  It  HE 
M=  1 

2C1  CO  NTIMt.lt 

KK=1 

CALL  V‘tULFF(E<n2, Xu, MM, *M, ,IE=) 

03  it?  1=1, MN 
113  X3  (I)  =  Xl(3  ) 

M=  M+l 

IF  ( (M*  DT ) #  EG. ..  .5)  THEN 
A  =  A  +  C .5 
DO  2.2  K= 1 , 2 

a  a 

IFCOEC.EG..)  THEN 

DO  2  l  I=i,  NC 

J=IC(I)  ' 

AA=  U  +  MODE  (R,  J)  *•  XI  (I) 

DO  2  »  1=1, NS 
J=IS  (I) 

A  A  =  aA+  VQ3E(K,J)  -*  X 1  ( I  ♦  n  C*  L  ) 

DO  2  f  1=1, NK 
J= In (I ) 

AA=  AA+  t'ODE  (  K  ,  J)  ^  X 1  ( I  +  nC’  A  +  NS  f  2 ) 

EL  SE 

no  2 5.1  1=1, NC 
J=  TC  < I) 

AA  =  AA  ♦  vcLE  «,J)  *  XKD 
00  212  1=1, NS 
j=ism 

AA  =  A 6  ♦  R  C  0  E  ( K  ,  J )  *  XI  (I  +  NC*.> 

E  M  01  r 
T(N)-  A  A 
f  J=  N  ♦  1 
CO  AT  It  UE 
IE  (A  «  Gr  #  2  ■  <i  )  GOTO  21* 

GOTO  2  3 
EL  SE 

GOTO  2'  1  .•) 

ENOI- 

21“  CONTINUE  ’* 

A  =  «> 

L=2M2  /A) 

PRINT*, *  TIME  X 

DO  1  1=1, L, Z 
W'TTEr  TU)  ,  T  ( I  +  i ) 


FORMAT  (  2X ,  F«  ,1  ,  <-X,FiL»6,-X,Fli>»  j ) 
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2  (. 

2  l 

Z.i 

211 

212 

2  2 


SUX: 


•■*■»** 


•J5F.0U7  INE  FCRPA  7k/ 7k  CPT  =  2 

l 

I 

I  SU EROUT INI  FGRvA{A,0,W,n,NZ,IC) 

CG  MM  Of./  M-  j.  N  1  /  f;  C  i  M 
REAL  /(KriKjNEIR)  ,W(ic)  ,  D  (:<-) 

;  INTEGER  IC(K) , I, J,N, M 

00  i  1=1, N? 

00  1  J=  1 1  NO 
A  (I,  ))  =  ■  . 

;  1  CONTINUE 

00  2  1=1,  M 
N  =  I C ( I ) 

A  (  (I  +  r  ))  =  0<M) 

\  A  (I,  (1*N>>  =  I.; 

’  A  (  (I+N>  ,1)  =  -  (W(M)*-*  Z) 

2  CO  NT  It  OE 

RETURN 
ENO 


COROUTINE  FO?MC 


7*-/ 71* 


0  P  T  = ; 


SUBROUTINE  FORN.C  (C,PHI$,  ,’J  ,1,2  ,  NS  E  N  ,  I : 
CO  WM  Cl/t  t  INI /NCI*: 

REAL  C  (  N  D I  ,  r.‘L'H)  ,  PHIS(  N01M,  i.L1 1) 

I N  TEG  rR  iC(l  )  )  NS  EN  ,NjN2,i,J,‘:> 

00  1  2=1, NEE N 
00  1  J=1,N2 
C(  I,  J)  =  '. .  ■ 

CO  NT  j.'  MUE 
DO  2  I=1,NSEN 
DO  2  J=1,N 
U  = IC  ( J) 

oc,  J)=  fh:s(m,d 

CO  NT  2  NUE 

ENC 


,J5hOU7JML  FLRP9 


7  k/ 7L 


C  PT=  -i 


SUBROUTINE  FC'\P3  (  3,?hI, NO,  CC» 


co  mm or/f 


il/N 


REAL  r<NDIf  , nOIM)  ,PRI(t  j  I  N,  N  El!' ) 

INTEGER  I C  <  f  )  ,  f,4C7,N,K,i  ,  J,  IF 

DO  1  1=1, NO 

00  1  J= 1 , NAf  7 

3(1,  J)  =  .  .1 

CONTEND' 

00  2  1=1, M 
00  2  J=  1 ,  N £  f  T 
i-  ICC) 

■*  ( c  *: ) ,  j)  -  - 1  *  ( i,  j) 

CON V  'MI- 
END 


•?» 


, 
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J6R0UTINS  FCRNX 


74/7 


rTN  3.1+ 


i 

< 


i 


1 


2 


3 


i 

j 

f 

i 

t 


1 


2 


3 


gpt=: 


SURr.OUTif  E  F O'  NX  i  (XO,INIT) 

COMMON/NU-'l/lCt  i2)  ,  IS  (i2>  ,  IF  ( 12)  ,M0,  '1S,NR 
REAL  XO  ( •  . )  » IN  I T  (.,12) 

INI  ESR  ► 

00  1  1=1, NO 
•■1=  IC  (I) 

XO (I )  =  INIT (1  ,M) 

XO<I*NC>  =  IrvIT  (  2 ,  R ) 

XO  (1+2*  NO)  =  INIT  (3,ri) 

XO  (I  *7'  NO  =  1  NIT  M) 

DO  2  7=1, NS 
H= IS  (3 ) 

XO  (1+ONO)  =  1  MT(1,M) 

XO  (I  NO+NS)  =  INIT  (  2,  ‘'.J 


XO  (I 

+<  - 

NC 

+  2k 

NS ) 

=  In  IKS, 

K) 

XO  (I 

+  / » 

NO 

♦  3* 

NS ) 

=  INIT  ( •* , 

F) 

no  3 

1  = 

1, 

NR 

M=  IF 

(I) 

XO  (I 

*l  - 

NO 

♦  :  » 

NS) 

=  I.MIKi, 

>-} 

XO  (1+M*:  +NS*J.+Niv)  =  INIT  <2,  N) 

JEN  .0 _ _ 

SUBROUTINE  FORt-.X'j  (XO  ,  If  l  T) 

CO  IT  /  NUN/ 1 C  (  1  2)  ,  IS  (  12)  ,  1'  (i2),N0,NS,N\ 
REAL  XO  (‘-  .  )  ,  IN  I T  (•>  ,  1 2  ) 

INTEOF-  F 

no  i  i=i, no 

H=  10(1) 

XO (I )  =  INI* (1 ,M) 

XO(I+NO)  =  1M7(2,K) 

XOd+rC‘2)  =  INITIO,1') 

XO  (I +NC*  0 )  =  IMTdijM) 

00.  2  1  =  1, NS 
■1=  IS  (I) 

XO  (I  +  N0*  -  )  =  I N I T  ( 1 , M ) 

XO  (I+tC?  -  +  MS  )  =  I  J1T  (  2,K 
DO  3  1  =  1,  N: 

M=  IF (I) 

X0(l+f0‘-  +  N  S  ~  2  )  =INIT(1,M) 

XO  (I  +  T  0‘  2  +NR )  =  INIT  (?,<•) 

EM  ^ 


SUKOUTir  L  pf  N'.  (VA1  ,  N,M) 

COWMEN  /  f'A  IM/N(  P 

c.E al 

I  NT'  V  R  \  ,1,  J 
P\I  NT','  * 

00  1  7=1,’ 

F  K I  ‘ 1 7  *  ( ;  >  F  1  2.t)  * ,  (KIT  (3  ,  J) ,  J  =  l ,  M) 
0  0  NT  I  m  E 

I NT  •{///)  • 
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